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ABSTRACT 


A  discrete-event  simulation  of  a  stochastic  process,  a 
machine-repairman  model,  has  been  programmed  on  the  IBM 
Personal  Computer.  The  model  consists  of  three  helicop¬ 
ters,  of  which  two  are  in  service  and  one  is  in  cold 
standby,  with  an  option  of  one  or  two  repairmen. 

The  program  output  is  a  graphics  display  containing  a 
system  state-versus -time  graph,  a  table  of  statistics,  and 
animated  figures  that  illustrate  the  current  state  of  the 
system.  The  program  user  can  directly  observe  the  dynamics 
of  the  model  as  the  fixed-increment,  simulation  clock, 
advances.  The  user  has  the  option  of  changing  the  follow¬ 
ing  model  parameters:  helicopter  failure  rate,  repairman 
service  rate,  and  the  number  of  repairmen  to  employ. 
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INTRODUCTION 


I . 

During  recent  years,  the  versatile  microcomputer  has 
found  a  niche  in  many  of  the  professional  discliplines 
including  education.  For  difficult  and  complex  subjects, 
educators  are  exploring  the  possibilities  of  the  microcom¬ 
puter  as  a  teaching  aid  in  addition  to  its  more  traditional 
role  as  a  computing  machine.  Some  universities  now  use 
microcomputer  graphics  to  illustrate  principles  that  cannot 
be  fully  explained  in  a  lecture.  [Ref.  1] 

A  widely  used  application  of  computers  is  to  imitate  or 
simulate  stochastic  processes.  This  technique  usually 
involves  only  the  numerical  evaluation  over  time  of  the 
corresponding  model  with  the  intent  of  gaining  some  under¬ 
standing  of  the  process.  But  using  graphics  as  part  of  the 
output  in  a  simulation  provides  the  user  another 
perspective.  For  example,  observing  a  working  model--such 
as  a  birth-death  process--can  help  in  grasping  the  funda¬ 
mentals  of  the  model  and  provide  some  insights  into  the 
process  not  easily  seen  on  the  classroom  chalkboard. 

This  paper  describes  the  development  of  a  discrete- 
event,  graphic  simulation  of  a  specific  example  of  the 
machine-repairman  model.  The  simulation  gives  the  observer 
the  impression  of  real-time  passage.  A  graph  on  the  visual 
display  shows  the  history  of  the  system's  transitions  from 


state  to  state,  and  animated  figures  illustrate  the 
dynamics  of  the  model.  The  simulation  also  exhibits  a 
table  of  the  theoretical,  and  the  current,  estimated  values 
of  quantities  commonly  used  in  evaluating  queueing  systems. 

Section  II  of  this  paper  is  a  discussion  of  the 
machine-repairman  model,  the  assumptions  upon  which  the 
model  is  based,  and  the  methods  used  for  obtaining  the 
theoretical  and  estimated  values  of  the  quantities  of 
interest.  Section  III  covers  the  simulation's  implementa¬ 
tion  on  the  IBM  Personal  Computer  (IBM  PC).  Appendix  A  is 
a  user's  guide  which  contains  summary  inf-  about  the 
machine-repairman  model  and  instructions  for  operating  the 
program.  The  program  flowchart  and  the  program  listings 
are  found  in  Appendixes  B,  C,  D,  and  E. 
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1 1  .  THE  MACHINE-REPAIRMAN  MODELS 


A.  DISCUSSION 

The  simulation  described  in  this  paper  actually  con¬ 
tains  two  machine-repairman  models.  In  both  models  there 
are  three  machines,  represented  by  helicopters.  The  models 
are  distinguishable  by  the  number  of  repairmen;  there  is 
one  repairman  in  the  first  and  two  repairmen  in  the  second. 

The  nominal  mission  is  to  keep  two  helicopters  flying 
at  all  times.  The  third  helicopter  is  placed  in  cold 
standby  as  a  backup  to  the  other  two.  If  available,  a 
repairman  is  assigned  immediately  to  a  helicopter  upon  its 
failure.  If  a  repairman  is  not  free,  then  the  helicopter 
joins  a  queue  to  await  repair. 

Several  assumptions  are  made  to  simplify  the  models.  In 
the  two-repairmen  model,  repairmen  are  assumed  to  be 
equally  competent  and  to  taxe  the  same  mean  time  to  repair 
a  helicopter.  Both  models  assume  the  repair  time  to  be 
continuous;  that  is,  repair  parts  are  always  available,  so 


there 

is  no  delay  in 

repairing 

a  machine 

other  than 

the 

time 

required  by  the 

repairman 

to  perform 

the  work. 

The 

repairmen  work  independently  and  do  not  assist  one  another. 
Similarly,  the  models  assume  the  helicopters  operate  inde¬ 
pendently  and  have  the  same  service  life  with  the  same  mean 
time  to  failure.  By  assumption,  a  helicopter  is  not 
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language  subroutines,  generating  the  graphics  display,  and 
performing  the  tasks  selected  from  the  program  menu.  Within 
the  main  program  are  routines  for  computing  the  limiting 
values  of  the  tabulated  quantities  and  for  printing  them  on 
the  display  screen.  The  main  program  also  calculates  the 
initial  failure  times  for  the  first  entry  into  the  clock 
module  . 

The  clock  module  determines  the  next-event  time, 
whether  the  event  is  a  failure  or  the  completion  of  repair, 
and  which  helicopter  it  affects.  The  simulation  clock  is  a 
fixed-increment,  time-advance  type.  In  contrast  to  a  next- 
event,  time-advance  clock  which  skips  over  the  intervening 
time  between  events,  the  clock  module  in  MACHREPR  advances 
the  time  in  equal  increments.  With  this  type  of  clock,  the 
user  experiences  the  illusion  of  real-time  passage. 

During  the  execution  of  the  simulation,  the  clock 
module  monitors  the  keyboard  for  certain  commands  issued  by 
the  user.  Pressing  the  "P"  key  pauses  the  program  to  allow 
the  user  to  study  the  graphic  display.  Pressing  the  "C" 
key  causes  the  program  to  continue.  The  "S"  key  will  stop 
the  simulation  and  return  the  program  to  the  program  menu. 
Should  the  clock  reach  the  preset  limit  of  9950  time  units, 
the  program  will  automatically  stop  the  simulation  and 
return  to  the  program  menu. 

Every  time  the  clock  is  advanced  one  increment,  the 
clock  module  updates  the  system  state-versus-time  graph 
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seen  above  the  others  is  in  active  service;  the  two 
machines  with  broken  rotors  are  down.  In  this  example,  the 
helicopter  graphics  indicate  that  the  system  is  in  state  2, 
where  helicopter  B  is  flying,  helicopter  C  is  under  repair, 
and  helicopter  A  is  awaiting  repair.  The  graph  shows  that 
the  current  time  is  69  and  that  the  system  jumped  from 
state  1  to  2  at  t  =  68. 

The  table  in  the  lower  right  corner  contains  the  esti¬ 
mated  and  the  limiting  values  of  the  quantities  discussed 
in  Section  II.  The  row  labels  are: 

Ar  -  the  arrival  rate  of  machines  into  the  repair 

queue , 

Wq  -  the  average  wait  time  of  a  machine  in  the  repair 
queue , 

Lq  -  the  average  number  of  machines  in  the  repair 

queue , 

W  -  the  average  amount  of  time  a  machine  is  down, 

L  -  the  average  number  of  machines  in  the  repair 

system, 

Avl  -  the  proportion  of  time  that  at  least  one  machine 
is  available. 

Pd  -  the  probability  that  a  down  machine  must  wait  for 
repair , 

Ps  -  the  probability  that  an  up  machine  must  go  to 
cold  standby. 

MACHREPR  consists  of  four  major  parts:  (1)  main  pro¬ 
gram,  (2)  clock  module,  (3)  failure  module,  and  (4)  repair 
module.  The  main  program  handles  the  housekeeping  details 
such  as  initializing  variables,  loading  the  machine 
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generator  seed,  or  (5)  end  the  program.  The  instructions 
that  are  available  on  the  screen  are  a  summarized  version 
of  those  found  in  Appendix  A,  the  user's  guide.  The  model 
parameters  that  can  be  changed  by  the  user  are  the  mean 
time  to  failure,  the  mean  time  for  repair,  and  the  number 
of  repairmen. 


Figure  3:  An  Example  of  the  Graphics  Display. 


Figure  3  shows  an  example  of  the  graphics  display 
screen  as  it  appears  during  the  simulation.  The  upper  half 
of  the  screen  contains  the  system  state-versus -time  graph. 
The  system  state  values  are  on  the  ordinate  of  the  graph, 
and  fifty  divisions  representing  time  are  on  the  abscissa. 
The  three  graphic  figures  of  helicopters  in  the  lower  left 
corner  symbolize  the  status  of  the  system.  The  helicopter 


B.  THE  PROGRAMMING  LANGUAGES 


MACHREPR  is  written  in  Advanced  BASIC  which  is  the  pro¬ 
gramming  language  supplied  with  PC-DOS.  Advanced  BASIC  was 
selected  because  of  its  convenient  graphics  commands  and 
its  ready  availability.  Although  Advanced  BASIC  is  an 
interpreter,  program  speed  was  not  an  essential  criterion 
that  warranted  using  a  compiled  language.  Advanced  BASIC 
was  insufficient,  however,  in  two  aspects,  and  so  MACHREPR 
uses  two  subroutines  written  in  8088  assembly  language  to 
improve  program  performance.  These  two  subroutines,  named 
SCRNSHFT  and  RNGEM ,  will  be  described  later  in  this 
section . 

C.  PROGRAM  DESCRIPTION 

Appendix  B  contains  the  flowcharts  for  the  BASIC  pro¬ 
gram  MACHREPR.  The  flowcharts  illustrate  the  simulation's 
general  construction  and  program  flow;  however,  they  are 
not  a  verbatim  duplication  of  the  program  listing.  The 
program  listing  for  MACHREPR  and  the  two  source-code  list¬ 
ings  for  RNGEN  and  SCRNSHFT  are  in  Appendixes  C,  D,  and  E 
respectively . 

After  initializing  several  variables  and  arrays,  the 
program  presents  the  title  screen.  The  program  menu 
appears  next  with  the  following  options:  (1)  see  the  pro¬ 
gram  instructions,  (2)  change  the  model's  parameters,  (3) 
use  the  default  parameters,  (4)  set  the  random-number 


III.  THE  SIMULATION 


A.  THE  COMPUTER 

The  microcomputer  used  for  this  simulation  was  the  IBM 
Personal  Computer.  Since  the  program  MACHREPR  uses  color 
graphics  extensively,  the  computer  requires  a  color  display 
monitor  and  either  IBM's  Color/Graphics  Adapter  or  a  suit¬ 
able  graphics  card  from  another  manufacturer.  The  program 
will  not  run  without  a  graphics  adapter.  Other  minimum, 
hardware  requirements  for  this  simulation  are  128  kilobytes 
of  read-write  memory  (RAM)  and  at  least  one  5i-inch,  floppy 
disk-drive . 

At  the  beginning  of  this  project,  the  author  chose  to 
use  IBM's  disk  operating  system  PC-DOS  with  the  intent  of 
enabling  the  simulation  to  run  on  other  microcomputers 
using  the  more  generic  version  of  PC-DOS:  MS-DOS  by 
Microsoft.  It  soon  became  apparent,  however,  that  MACHREPR 
would  need  to  directly  access  certain  parts  of  the 
computer's  memory  in  order  to  improve  the  program's 
performance.  Since  the  program  makes  direct  calls  to  the 
PC-DOS,  there  can  be  no  assurance  that  the  simulation  will 
"IBM-compatible"  machines. 


run  on  other 


For  P  ,  the  only  time  a  machine  must  wait  in  cold  standby 
is  when  the  system  is  in  state  0;  therefore,  Ps  =  p^  . 

The  availibility  of  at  least  one  helicopter  in  the 
system  is  the  sum  over  all  states  n  of  the  percentage  of 
the  flyable  helicopters  in  each  state  multiplied  by  the 
proportion  of  time  the  system  is  in  that  state; 

Av  =  p  q ( 1 )  +  p^(2/3)  +  p2(l/3)  +  p3(0). 

The  estimate  of  the  availabiliy  is  calculated  by  substitut¬ 
ing  p  ,  n  =  0,1, 2, 3,  into  the  above  expression. 


r- ; 


w  v  -  v  ■  -v  -  v  ■  "  'j  ' 


rt’VVV  *.■  •  ’ 


B 


The  average  waiting  time  quantities,  W  and  are  then 
calculated  by  using  Little's  formulas 


W  =  L/X 
"q  = 


To  compute  the  estimate  W  ,  the  program  keeps  a  record  of 
the  length  of  time  each  helicopter  spends  in  the  repair 
queue  and  uses  the  expression 


W  =  Cl/m]  J  q. 
q  ikl  1 


m 


m  =  0,1,2,.. 


where  q^  is  the  waiting  time  in  queue  of  the  ith  helicopter 
and  m  is  the  total  number  of  helicopters  that  have  entered 
the  queue  by  time  t.  Then  W  is  easily  found  with 

W  =  Wq  +  E ( S  ) . 

The  estimate  of  the  probability  that  a  down  machine 
must  wait  for  repair  is  the  proportion  of  the  time  that  the 
system  spends  in  a  state  in  which  a  machine  is  down  with  no 
unemployed  repairman.  The  repairman  is  not  available  in 
the  one-repairman  model  when  the  system  is  in  state  2  or  3: 

Pd  =  p2  +  Pa¬ 
in  the  two-repairmen  case,  the  repairmen  are  busy  when  the 
system  is  in  state  3: 

Pd  =  P3- 
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In  the  two-repairmen  case,  a  queue  forms  only  when  the 
system  is  in  state  3,  thus 


I 

i 

I 


< 


(• 


j 


i 


and 


The  average  rate  at  which  helicopters  enter  the  repair 
queue,  X,  is  dependent  upon  the  number  of  active  machines 
subject  to  failure  in  state  n  and  the  proportion  of  time 
the  system  is  in  state  n.  There  are  two  active  helicopters 
in  states  0  and  1,  one  in  state  2,  and  none  in  state  3.  The 
combined  failure  rate  of  the  active  helicopters  in  state  n 
is  given  by 


a 


n 


2a  n  =  0,1 

<  a  n  =  2 

0  n  =  3 . 


Then  X  is  found  by 


X 


n^O 


U  anpn 


2ap0  +  2 ap1  +  ap2 . 


For  the  estimate  "x ,  the  program  divides  the  number  of  fail¬ 
ures  that  have  occurred  in  time  t  by  t,  that  is 

I 


X  =  total  number  of  failures/total  elapsed  time. 


failure  rate  and  the  repairman  service  rate,  will  not  be 
discussed  in  this  paper.  The  interested  reader  should 
refer  to  the  bibliography  for  a  list  of  textbooks  on  queue¬ 
ing  theory  and  birth-death  processes.  In  addition  to  the 
limiting  probabilities,  {p  ,  n  =  0,1, 2, 3},  the  program 
computes  the  fraction  of  the  time  the  system  is  in  state  n, 
denoted  p  ,  as  equal  to  the  total  time  in  state  n  divided 
by  the  total  elapsed  time.  The  quantity  pn  is  an  estimate 
of  pn. 

The  value  of  L,  the  theoretical  average  number  of 
machines  in  the  repair  system,  is  found  by  averaging  the 
number  of  down  machines  over  all  states: 


.3 


L  =  )  np  =  (l)p,  +  (2)p?  +  (3)p 

nio  J 


The  estimate  of  L,  denoted  L,  is  computed  with  the  same 


expression  with  the  exception  that  the  estimate  is  sub¬ 


stituted  for  p 


The  average  number  of  machines  in  the  repair  queue  L 


and  its  estimate  depends  in  part  upon  the  number  of 


repairmen  in  the  model.  In  the  one-repairman  case, 
machines  occupy  the  repair  queue  only  when  the  system  is  in 


state  2  or  3 ,  thus 


As  shown  m  Figures  1  and  2,  the  variable  q  is  the  time 
that  a  down  machine  waits  in  the  repair  queue.  Then 
W  =  E(q).  It  follows  that  the  average  total  time  that  a 

q 

helicopter  is  in  the  repair  system  is  the  sum  of  the 
average  waiting  time  in  queue  and  the  expected  service 
time:W=W  +E(S). 

q 

Little's  formulas  show  the  relationship  of  the  four 
quantities  of  interest 


L  =  XW 

q  q 

where  X  is  defined  as  the  average  arrival  rate  of  machines 
entering  the  repair-queue  system. 

There  are  three  additional  quantities  that  will  be 
computed  and  presented  by  the  simulation.  They  are  defined 


=  P{down  machine  must  wait  for  repair}, 

=  P{up  machine  must  go  to  cold  standby}, 

=  availability  -  the  proportion  of  time  that  at 
least  one  helicopter  is  not  down. 


To  calculate  all  of  eight  of  these  values,  it  is  first 
necessary  to  compute  the  system's  steady-state  probabili¬ 
ties,  p^  =  P  { N  =  n},  n  =  0,1,2, 3.  In  other  words,  pn  is 
the  long-run  probability  that  the  system  is  in  state  n.  The 
details  of  calculating  p^ ,  a  function  of  the  machine 


of  a  helicopter  until  failure  and  assume  that  T  is  exoonen- 
tiallv  distributed  with  a  failure  rate  a.  Then  the  mean 
time  to  failure  (MTTF)  of  a  helicopter  is  E(T)  =  1/a. 
Similarly,  the  completion  of  repair  on  a  machine  is  a 
death,  so  define  S  to  be  the  service  time  of  a  repairman. 


The  service 

time  is 

also  assumed  to 

have 

an  exponential 

distribution 

with  a 

service  rate  p. 

The 

mean  time  for 

repair  (MTFRj  is  then  E(S)  =  l/p. 

For  these  machine-repairman  models,  the  state  of  the 
system,  N(t)  =  n,  is  defined  as  the  number  of  helicopters 
down  at  time  t.  Since  there  are  only  three  helicopters  in 
the  models,  then  n  =  0,1, 2, 3. 


B.  MODEL  ESTIMATORS 

In  the  study  of  queueing  models,  there  are  four  quanti¬ 
ties  that  are  commonly  used  to  evaluate  the  queueing  sys¬ 
tem's  performance.  As  part  of  the  graphical  output,  the 
simulation  will  provide  the  user  a  table  of  theoretical  and 
estimated  values  of  these  four  quantities: 


w 

w 


the  average 
system. 

number 

of 

machines 

in 

the 

repair 

the 

average 

number 

of 

machines 

in 

the 

repair 

queue , 

the  average  amount  of  time  a  machine  is  down, 

the  average  amount  of  time  a  machine  waits  for 
repair . 
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The  cyclic  nature  of  the  model  is  clearly  evident  in 
the  figure.  A  helicopter's  life  cycle  consists  of  an  oper¬ 
ating  period  followed  by,  after  a  possible  delay  in  queue, 
a  repair  period  which  is  in  turn  followed  by,  after  another 
possible  delay  in  cold  standby,  another  operating  period. 


Figure  2:  Machine-Repairman  Model  with  Two  Repairmen. 

Figure  2  depicts  a  pictorial  representation  of  the 
machine-repairman  model  with  two  repairmen  [Ref.  3].  In 
this  system,  helicopters  waiting  in  the  repair  queue  go  to 
the  first  available  repairman  for  service.  The  maximum 
possible  length  of  the  repair  queue  is  only  one.  As  in  the 
case  of  the  one-repairman  model,  the  standby  queue  can 
contain  at  most  one  helicopter. 


The  machine-repairman  model  is  a  special  case  of  a 
birth-death  process  in  which  the  event  of  a  machine  failing 
at  time  t  constitutes  a  birth.  Let  T  be  the  operating  life 


that  is  down  and  awaiting  repair  is  colored  red.  When  the 
down  machine  is  in  repair,  it  is  seen  with  a  repairman. 

Figure  1  is  a  frequently  used,  pictorial  representation 
(see  [Ref.  2],  for  instance)  of  the  machine-repairman  model 
with  one  repairman.  The  dashed  lines  enclose  the  model's 
repair-queue  system  which  contains  a  single-server  queue. 
The  machines  outside  the  dashed  box  are  the  potential 
customers  for  the  repair  queue.  Since  there  are  only  three 
machines  in  the  model,  the  longest  possible  length  of  the 
queue  will  be  two. 


Standby 


Helo 


Helo 


Repair  |  Repairman 


Figure  1:  Machine-Repairman  Model  with  One  Repairman, 


This  model  has  a  second  queue,  the  cold-standby  queue, 
since  a  helicopter  may  have  to  wait  before  being  assigned 
to  a  mission.  The  standby  queue  can  contain  only  one  heli¬ 
copter,  and  this  occurs  when  there  are  no  machines  in  the 
repair-queue  system. 


subject  to  failure  while  in  cold  standby  since  it  is  not 
operating . 

A  machine-repairman  model  typically  assumes  its 
machines  operate  continuously  until  failure.  In  the  case 
of  aircraft,  this  assumption  does  not  reflect  reality  since 
in  actual  practice  aircraft  do  not  fly  continuously,  but 
have  instead  periods  of  inactivity  between  missions  of 
limited  duration.  As  already  noted,  one  can  reasonably 
assume  that  an  inactive  helicopter  is  not  subject  to  fail¬ 
ure  and  that  the  failure  process  occurs  only  when  the 
helicopter  is  active. 

The  actual  failure  rate  of  a  flying  helicopter  is  the 
number  of  failures  per  flight  hour.  The  product  of  the 
actual  failure  rate  times  the  ratio  of  the  number  of  flight 
hours  performed  to  the  number  of  calendar  hours  elapsed 
until  failure  defines  the  helicopter's  effective  failure 
rate.  Thus  the  effective  failure  rate  is  the  number  of 
failures  per  calendar  hour.  To  simplify  the  simulation, 
the  models  assume  that  the  helicopters  operate  continuously 
until  failure  with  the  effective  failure  rate. 

The  selection  of  helicopters  to  represent  the  machines 
in  these  models  allows  a  simple  graphic  design  in  depicting 
their  status.  An  operating  machine  is  seen  in  the  display 
as  a  green,  flying  helicopter;  a  helicopter  that  is  in  cold 
standby  is  also  colored  green  but  is  not  flying.  A  machine 
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according  to  the  current  state  and  time  values.  As  the 
time  advances,  the  plotted  lines  will  eventually  reach  the 
edge  of  the  graph.  When  this  condition  is  met,  the  clock 
module  shifts  the  entire  graph  left  one  time  unit  and  draws 
the  new  time  increment.  This  procedure  continues  for  all 
successive  time  increments.  When  the  simulation  clock  has 
reached  the'  time  for  the  next  event,  the  program  jumps  to 
either  the  failure  module  or  the  repair  module  as 
appropriate . 

The  failure  module  changes  the  helicopter  display 
according  to  the  current  state  and  model  parameters.  The 
program  keeps  track  of  the  machines  individually,  so  that 
the  graphic  image  of  the  particular,  failing  helicopter 
transforms  from  a  green,  flying  figure  to  a  red,  grounded 
one.  If  a  repairman  is  free,  the  graphic  figure  of  a  man 
is  drawn  next  to  the  aircraft.  If  no  repairman  is  avail¬ 
able,  then  the  human  figure  is  omitted,  signifying  that  the 
helicopter  is  placed  in  the  repair  queue.  If  a  repairman 
was  assigned,  the  failure  module  calls  the  random-number 
generator  and  computes  the  service  time.  The  sum  of  the 
service  time  and  the  current  time  is  the  time  that  the 
repair  will  be  completed  and  this  is  compared  to  the  next- 
event  times  of  the  other  machines  after  the  program  returns 
to  the  clock  module. 

In  a  similiar  manner,  the  repair  module  changes  the 
graphic  figure  of  the  helicopter  that  has  just  completed 
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repairs.  If  there  are  less  than  two  helicopters  flying 

when  the  repair  event  occurs ,  then  the  repaired  machine 
transforms  from  red  colored  and  disabled  to  green  colored 
and  flying.  If  two  helicopters  are  already  flying,  then 
the  repaired  machine  goes  into  the  standby  queue.  The  now 
free  repairman  moves  to  t'he  helicopter  at  the  head  of  the 
queue,  or  if  the  queue  is  empty,  he  disappears  from  the 

screen . 

Before  the  program  returns  to  the  clock  module,  the 
repair  module  computes  the  interarrival  time  to  failure  of 
the  helicopter  just  repaired.  Failure  times  are  not  com¬ 
puted  for  helicopters  placed  in  the  standby  queue.  If  the 
repairman  is  assigned  to  another  helicopter,  the  repair 

module  also  computes  the  service  time. 

After  returning  to  the  clock  module,  the  simulation 

updates  the  estimates  of  the  tabulated  quantities  and 
prints  them  on  the  display  screen.  Both  the  failure  and 
repair  modules  change  the  state-indicator  variable  to  its 
new  value  so  that  the  clock  module  will  draw  a  horizontal 
line  one  time  unit  in  length  at  the  proper  place  on  the 
graph . 

SCRNSHFT  is  the  subroutine  that  shifts  the  system 
state-versus -time  graph  one  time  unit  to  the  left  per  clock 
cycle.  The  subroutine  uses  certain  program  transfers 
called  interrupts  to  the  PC-DOS  Basic  Input/Output  System 
(BIOS)  which  contains  routines  that  control  the  video 
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display.  SCRNSHFT  is  not  appreciably  faster  than  a  BASIC 
routine  designed  to  perform  the  same  task;  however,  the 
subroutine  eliminates  an  annoying  flicker  in  the  display 
that  was  characteristic  of  the  BASIC  version. 

D.  PROGRAMMING  CONSIDERATIONS 

The  program  computes  the  interarrival  (failure)  and 
service  times,  assumed  to  be  continuous,  exponentially  dis¬ 
tributed,  random  variables,  in  the  single-precision  format. 
The  current  time,  an  integer  value,  is  added  to  the  inter¬ 
arrival  time  and  to  the  service  time  to  yield  the  next- 
failure  time  and  the  next-repair  time  respectively.  In  the 
unlikely  event  that  the  two  single-precision,  random 
variables  match,  the  simulation  will  proceed  as  if  the 
failure  occurred  earlier. 

Whenever  the  model  jumps  from  one  state  to  another,  the 
program  computes  a  next-failure  time  and  a  next-repair  time 
for  the  helicopter  just  affected  by  the  transition.  Thus 
MACHREPR  always  keeps  in  memory  six  next-event  times  for 
all  machines  regardless  of  their  individual  status.  To 
determine  the  next  event,  the  clock  module  compares  the 
smallest  of  the  three  failure  times  with  the  smallest  of 
the  three  repair  times.  With  this  method,  there  must  be  a 
means  of  preventing  the  program  from  simulating  events  con¬ 
sidered  impossible  by  the  assumptions  of  the  model. 
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For  example,  one  model  assumption  is  that  a  helicopter 

in  cold  standby  is  not  subject  to  failure.  It  follows  that 

the  same  machine  cannot  be  repaired  since  it  is  not  in  the 

repair  queue.  Similarly,  a  helicopter  that  is  in  repair 

cannot  fail  and  a  helicopter  that  is  flying  cannot  be 

repaired.  To  handle  these  impossible  events,  the  program 

will  set  the  failure  time  or  the  repair  time,  as  appropri- 
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ate,  to  equal  10  which  represents  infinity. 

Because  of  the  discrete  nature  of  the  graphic  display, 
all  time  values  must  be  converted  to  integers  prior  to 
plotting  the  state-versus-time  graph.  It  is  for  this  rea¬ 
son  that  the  random  number  representing  the  next-event 
time,  be  it  a  failure  or  a  repair,  is  rounded  up  to  the 
next  higher  integer.  The  result  is  that  the  interarrival 
and  service  times  are  not  exponentially  distributed  but  are 
instead  geometrically  distributed. 

The  discreteness  of  the  time  variable  introduces  a  bias 
in  the  estimators  since  they  are  based  on  the  assumption  of 
continuous  distributions.  The  effects  of  the  bias,  how¬ 
ever,  are  minimal.  The  values  of  the  estimates  are  still 
seen  to  converge  to  the  limiting  values  as  time  advances. 
Despite  the  generation  of  geometric  random  variables,  the 
simulation  still  provides  a  sufficient  visual  demonstration 
of  the  machine-repairman  model. 
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E.  RANDOM  NUMBER  GENERATION 

The  author  found  BASIC's  random-number  generator,  a 
function  called  RND ,  unsuitable  for  simulation.  Even  so,  a 
sample  of  10,000  uniform  variates  produced  by  RND  passed 
with  .95  confidence  three  standard,  nonparametr ic  tests  for 
uniformity  and  independence:  the  serial  test,  the  frequency 
test,  and  the  runs  up-and-down  test.  On  the  other  hand,  a 
two-dimensional  scatter  plot  of  5000  pairs  of  uniform  vari¬ 
ates,  shown  in  Figure  4,  clearly  illustrate  the  lattice 
structure  of  the  RND  generator. 


Figure  4:  Two-Dimensional  Plot  of  5000  Pairs  of  Uniform 
Variates  from  the  BASIC  RND  Function. 


MACHREPR  uses  a  random-number  generator  called  RNGEN . 
Written  in  8088  assembly  language,  RNGEN  is  modification  of 


a  generator  programmed  by  Associate  Professor  B.O.  Shubert 
of  the  Naval  Postgraduate  School  for  use  on  the  IBM  PC. 
RNGEN  uses  the  algorithm  described  by  Fishman  [Ref.  4] 
which  is  based  on  the  expression 

Xi+1  =  16807*Xi  mod  (231  -  1) 

where  X.  is  the  old  seed  and  X.  is  the  new  seed. 
1  l+l 

Dividing  the  new  seed  by  the  modulus  produces  the  uniform 
(0,1),  random  number. 

This  algorithm  is  the  one  used  by  LLRANDOM,  a  random- 
number  generator  developed  at  the  Naval  Postgraduate  School 
and  described  by  Lewis,  Goodman,  and  Miller  [Ref.  5]  as 
able  to  produce  good  quality,  uniform  variates. 

The  simulation  computes  exponential  random  numbers  by 
the  inverse  probability  integral  transformation  of  uniform 
variates  produced  by  RNGEN.  Let  U  be  a  uniformly  distri¬ 
buted,  random  number  and  let  X  be  exponentially  distributed 
with  parameter  a.  An  exponential  random  number  is 
generated  by  setting  u  =  F(x)  and  solving  for  x, 

u  =  F  ( x )  =  1  -  e 

.  -  ax 

1  -  u  =  e 

ln(l-u)  =  -ax 

x  =  -ln(l-u)/a,  or 

x  =  -ln(u)/a. 
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A.  INTRODUCTION 


MACHREPR  is  a  stochastic,  discrete-event  simulation  of 
a  machine -repairman  model,  which  is  a  special  case  of 


birth-death 

processes . 

The 

model  consists 

of 

three  heli- 

copters ,  of 

which  two 

are 

in  service  and 

one 

is  in  cold 

standby,  with  an  option  of  one  or  two  repairmen.  The  pro¬ 
gram  output  is  a  graphics  display  composed  of  a  system 
state-versus -time  graph,  a  table  of  statistics,  and  animat¬ 
ed  figures  that  illustrate  the  dynamics  of  the  process. 

This  manual  discusses  the  model  and  guides  the  user  in 
the  step-by-step  operation  of  the  program.  In  describing 
commands  and  keyboard  entries,  the  following  notation  is 
used.  Single  quotes  ('  ')  enclose  commands  and  phrases 
that  are  to  be  entered  verbatim.  The  act  of  pressing  a 
specific  key  is  indicated  by  brackets  (<  >)  encasing  the 
letter  or  letters  on  the  key.  For  example,  'format  b:' 
<  ENTER  >  means  type  the  command  exactly  as  it  appears 
between  the  quotes  and  follow  it  by  pressing  the  enter  key. 
This  user's  guide  shows  the  commands  in  the  lower  case; 
however,  the  IBM  PC  will  accept  upper  case  letters  as  well. 

B.  MODEL  DESCRIPTION 

MACHREPR  contains  two  versions  of  the  machine-repairman 
model:  the  one -repai  rman  case  and  the  two-repairman  case. 
In  both  cases,  the  nominal  mission  of  the  model  is  to  keep 
two  helicopters  flying  at  all  times.  If  the  third 
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helicopter  is  serviceable,  it  is  kept  in  cold  standby  as  an 
immediate  replacement  for  the  other  two.  Upon  failure  of  a 
helicopter,  a  repairman  is  assigned  if  he  is  available; 
otherwise  the  helicopter  joins  a  queue  to  await  repair. 

Disabled  (down)  helicopters  remain  in  the  repair  queue 
until  a  repairman  becomes  free.  The  repair  queue  is  a 
first  in,  first  out  type  (FIFO).  A  newly  repaired  helicop¬ 
ter  will  go  immediately  into  flying  operation  if  less  than 
two  aircraft  are  in  service.  If  there  are  already  two 
helicopters  flying  when  a  repair  is  completed,  then  the 
third  goes  into  the  cold-standy  queue. 

There  are  several  simplifying  assumptions  for  these 
models.  All  helicopters  are  assumed  to  have  exponentially 
distributed  service  lives  with  the  same  failure  rate  a.  If 
T  is  defined  as  the  length  of  a  helicopter's  operating 
life,  then  the  mean  time  to  failure  ( MTTF )  is  E(T)  =  1/a. 
The  simulation  uses  the  MTTF  as  the  input  parameter  that 
establishes  the  particular  exponential  distribution  for  the 
operating  life.  Helicopters  will  fly  continuously  until 
failure.  A  machine  that  is  in  cold  standby  is  assumed  not 
to  be  subject  to  failure. 

Another  assumption  is  that  the  service  time  for  the 
repairmen  has  an  exponential  distribution  with  parameter  u. 
If  S  is  the  time  required  to  complete  a  repair  job,  then 
the  mean  time  for  repair  is  E(S)  =  l/u.  Like  the  MTTF,  the 
mean  time  for  repair  ( MTFR )  is  an  input  variable  for  the 
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simulation.  The  repairmen  are  assumed  to  be  equally  compe¬ 
tent  and  thus  have  the  same  MTFR. 

In  both  cases  of  the  machine-repairman  model,  the  state 
of  the  system  is  defined  as  the  number  of  helicopters  down 
at  time  t.  The  state  variable  is  N(t)  =  n,  n  =  0,1, 2, 3. 

C.  HARDWARE  AND  SOFTWARE 

1 .  Hardware 

The  simulation  is  programmed  for  the  IBM  Personal 
Computer  (IBM  PC)  with  a  color  monitor.  Since  the  program 
uses  color  graphics,  the  microcomputer  must  be  equipped 
with  IBM's  Color/Graphics  Adapter  or  a  suitable  graphics 
card  from  another  manufacturer.  MACHREPR  will  not  run 
without  the  graphics  adapter.  The  program  uses  two 
assembly-language  subroutines  that  are  loaded  in  memory 
outside  of  BASIC's  64  kilobyte  (K)  workspace;  therefore, 
there  must  be  a  minimum  of  128  K  of  memory  installed  in  the 
microcomputer.  The  computer  must  have  at  least  one 
5i-inch,  floppy  disk-drive. 

2 .  Software 

MACHREPR  and  the  two  subroutines  are  supplied  on  a 
distribution  diskette.  The  user  must  provide  the  disk 
operating  system  PC-DOS  2.0  (or  2.10)  and  Advanced  BASIC 
(BASICA)  in  order  to  run  the  simulation.  MACHREPR  uses 
BASIC's  POKE  command  and  certain  program  transfers,  called 
interrupts,  to  the  Basic  Input/Output  System  (BIOS)  in  the 


IBM  ROM  (read-only  memory).  Because  of  these  interrupts, 
there  can  be  no  assurance  that  the  simulation  will  operate 
properly  on  other  "IBM  compatable"  machines. 

3 .  Program  Files 

The  distribution  diskette  contains  the  following 
program  files: 

*  MACHREPR . BAS  -  the  main  program. 

*  RNGEN.SRT  -  a  random-number  generator  written  in  8088 
assembly-language  used  as  a  subroutine  to  the  main 
program. 

*  SCRNSHFT . SRT  -  another  assembly-language  subroutine 
used  to  perform  a  specific  task  on  the  graphics 
screen . 

*  RNGEN.LST  -  a  source-code  listing. 

*  SCRNSHFT. LST  -  a  source-code  listing. 

*  MACHREPR .  BAT  -  a  batch  file  that  loads  and  starts  the 
simulation . 

The  simulation  requires  the  first  three  files  listed  above 
in  order  to  run.  The  two  subroutines  are  binary-image  files 
and  can  not  be  executed  from  PC-DOS.  The  two  source-code 
listings  are  provided  as  a  convenience  to  the  user. 
MACHREPR .  BAT  is  a  simple  example  of  how  a  batch  file  can  be 
used  to  easily  start  the  simulation. 

D.  GETTING  STARTED 

1 .  Making  a  Backup  Copy 

The  first  step  is  to  make  a  copy  of  the  distribu¬ 
tion  diskette.  After  the  copy  is  made,  the  original 


38 


graph.  With  exponentially  distributed  interarrival  and 


service 

sojourn 


times,  there  is  a  non-zero  probability 
in  any  particular  state  will  be  less  than 


that 

one 


unit.  The  program  will  still  alter  the  helicopter  display 
properly  and  compute  the  estimates  correctly;  however,  the 


state-versus -time  graph  may  display  a  jump  of  two  states 
because  of  the  program's  inability  to  plot  a  line  of  less 
than  one  time  unit  in  length.  The  number  of  times  this 
discrepancy  occurs  can  be  minimized  if  the  MTTF  and  the 


MTFR  are  set  at  values  greater  than  ten. 
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a  result,  the  estimators  are  biased  since  they  are  based  on 
the  assumption  of  continuous  time  when  in  fact  they  are 
calculated  with  data  based  on  discrete  increments  of  time. 


The  effects  of  the  bias  are  minimal;  the  simulation  still 
provides  a  sufficient  demonstration  of  the  convergence  of 
the  estimates  to  the  limiting  values. 

G.  KEYBOARD  COMMANDS 

There  are  three  keyboard  commands  available  to  the  user 
as  the  simulation  is  in  progress.  When  <p>  (for  pause)  is 
pressed,  the  simulation  freezes  to  give  the  user  the  oppor¬ 
tunity  to  study  the  graphic  display.  To  continue  the  simu¬ 
lation,  press  <c> .  The  third  command  is  <s>  for  stop  the 
simulation.  The  <s>  key  will  send  the  program  back  to  the 
program  menu. 

H.  PROGRAM  LIMITATIONS 

MACHREPR  always  starts  a  simulation  under  the  same  con¬ 
ditions.  The  initial  state  is  N(0)  =  0  with  helicopters  B 
and  C  in  flight  and  helicopter  A  in  cold  standby.  There 
are  no  provisions  in  the  program  for  starting  the 
simulation  in  another  state.  Restarting  the  simulation 
from  the  program  menu  causes  the  current  time  and  all  data 
collection  variables  to  be  set  to  zero. 

As  noted  before,  all  event  times  are  rounded  up  to  the 
next  higher  integer  for  plotting  the  state-versus -time 
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*  Lq  -  The  average  number  of  machines  in  the  repair 

queue . 

*  W  -  The  average  amount  of  time  that  a  machine  is  in 

the  repair  system.  This  quantity  is  the  sum  of 
Wq  and  the  expected  service  time. 

*  L  -  The  average  number  of  machines  in  the  repair 

system.  This  quantity  includes  those  machines 

in  the  repair  queue  and  those  machines  in 

repair. 

*  Avl -  The'  machine  availability  defined  as  the 

proportion  of  the  time  that  at  least  one 
helicopter  is  flyable. 

*  Pd  -  The  probability  that  a  down  machine  must  wait 

for  repair. 

*  Ps  -  The  probability  that  an  up  machine  must  wait  in 

cold  standby. 

The  numbers  that  appear  in  the  right-hand  column 
under  the  heading  "Limit"  are  the  limiting  values  of  the 
quantities  of  interest.  In  other  words,  the  quantities 
converge  to  the  limiting  values  as  time  goes  to  infinity. 
The  limiting  values  are  calculated  with  the  balance  equa¬ 
tions  that  correspond  with  the  particular  machine-repairman 
model . 

The  values  in  the  left-hand  column  are  estimates  of 
the  limiting  values.  These  quantities  are  computed  with 
data  collected  from  the  simulation.  Upon  the  occurance  of 
an  event,  the  program  recomputes  the  estimates  after  taking 
into  account  the  newly  collected  data  from  the  last  state. 

Because  the  plots  on  the  state-versus-time  graph 
must  be  given  in  terms  of  integer  coordinates,  all  time 
variables  must  be  rounded  up  to  the  next  higher  integer.  As 
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Figure  11:  Helicopter  Display  in  State  3. 


3 .  The  Quantities  of  Interest 

In  the  study  of  queueing  systems,  there  are  several 
quantities  of  interest  that  are  used  to  evaluate  the  sys¬ 
tem's  performance.  In  the  lower  right  corner  of  the 
graphics  display  (see  Figure  7)  there  is  a  table  of  the 
more  commonly  used  quantities. 


The  row  labels  have  the  following  definitions: 


*  AR 

-  The  average 
repair  system 

arrival 

rate 

of  machines  into 

the 

*  Wq 

-  The  average 
repair  queue. 

waiting 

time 

of  a  machine  in 

the 

and  service  (repair)  times.  The  program  changes  the 
appropriate  graphic  figures  when  the  simulation  clock 
reaches  the  time  for  the  next  event. 

Figure  10  shows  an  example  of  a  helicopter  display 
that  indicates  state  2.  In  this  case,  C  is  in  repair. 
Helicopter  A  is  also  disabled,  and  since  it  does  not  have  a 
repairman,  A  is  in  the  repair  queue. 


Figure  10:  Helicopter  Display  in  State  2. 


In  Figure  11,  all  three  machines  are  down,  thus  the 
system  is  in  state  3.  There  is  a  repair  queue  of  length 
two  in  this  example. 

It  should  be  noted  that  Figures  10  and  11  are 
illustrations  from  the  one-repairman  model  since  only  one 
repairman  appears  in  either  display.  If  the  simulation  is 


Figure  8:  Initial  Helicopter  Display  in  State  0. 


Figure  9:  Helicopter  Display  in  State  1. 

As  the  simulation  proceeds,  random  numbers  are  gen¬ 
erated  and  transformed  into  interarrival  (failure)  times 
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at  the  height  corresponding  to  the  current  state  of  the 
system.  The  accumulation  of  these  individual  plots  forms  a 
step  graph  as  shown  by  the  example  in  Figure  7.  When  the 

graph  reaches  the  end  of  the  display,  the  program  makes 
room  for  the  next  plot  by  calling  SCRNSHFT  to  shift  the 

entire  graph  one  unit  to  the  left.  In  this  way,  the 

model's  activities  for  the  past  fifty  time  units  can  be 
seen  in  a  single  glance. 

2 .  The  Helicopter  Display 

The  three  helicopter  images  in  the  lower  left 
corner  of  the  graphic  display  indicate  the  current  state  of 
the  system.  The  display  in  Figure  8  shows  how  the  helicop¬ 
ters  are  arranged  in  the  initial  state  0.  All  three 

machines  are  colored  green.  Helicopters  B  and  C,  seen  in 

the  upper  half  of  the  display,  are  flying.  Helicopter  A  on 
the  "ground"  is  in  cold  standby. 

Since  all  helicopters  are  initially  either  flying 
or  in  cold  standby,  the  next  event  will  be  a  failure  in 

either  machine  B  or  C.  A  machine  that  fails  is  grounded 

and  transformed  into  a  red  colored  helicopter  with  a  broken 
rotor.  Because  there  is  only  one  down  helicopter,  the  sys¬ 
tem  is  in  state  1  and  so  there  is  at  least  one  idle 

repairman.  The  program  draws  a  graphic  figure  of  a  man 
beside  the  down  aircraft  to  represent  a  machine  that  is  in 
repair.  For  example.  Figure  9  depicts  state  1  with 


helicopter  B  in  repair. 


graphics  display  is  divided  into  three  sections:  the  state- 
versus-time  graph,  the  helicopter  display,  and  the  table  of 
quantities  used  with  queueing  models. 


Figure  7:  An  Example  of  the  Graphics  Display. 


1 .  The  State-Verses -Time  Graph 

The  upper  half  of  the  display  contains  the  state- 
versus-time  graph.  The  values  of  the  state  variable  N(t), 
the  number  of  machines  down  at  time  t,  are  on  the  ordinate. 
There  are  fifty  divisions  on  the  abscissa  which  represent 
units  of  time.  The  current  time  of  the  system  is  printed 
at  the  right-hand  end  of  the  abscissa. 

With  each  unit  advancement  of  the  simulation  clock, 
the  program  plots  a  horizontal  line  one  time  unit  in  length 


and  function  properly  with  non-integer  values  of  MTTF  and 
MTFR ;  however,  these  values  will  be  printed  on  the  graphics 
display  as  integers.  Negative  numbers  will  not  be  accepted 
for  the  MTTF  and  the  MTFR.  The  digits  1  and  2  are  the  only 
valid  responses  for  the  number  of  repairmen.  The  program 
will  immediately  start  the  simulation  after  the  number  of 
repairmen  has  been  entered. 

3  .  The  Default  Parameters 

If  this  option  is  selected,  then  the  simulation 
starts  with  the  model  parameters  set  at  default  values. 
These  values  are: 

MTTF  =  15  time  units, 

MTFR  =  10  time  units. 

Number  of  repairmen  =  1. 

The  default  parameters  are  selected  by  pressing  <d> . 

4 .  Ending  the  Program 

To  end  the  program,  press  <e>.  MACHREPR  will  clear 
the  screen,  exit  the  graphics  mode,  and  leave  the  computer 
in  the  BASIC  command  mode.  If  the  user  wants  to  exit  BASIC 
and  return  to  the  DOS,  type  'system'  <ENTER>. 


F.  THE  GRAPHICS  DISPLAY 

MACHREPR  is  designed  to  give  the  observer  the  illusion 
of  real-time  passage  as  the  probabilistic  events  occur  in 
the  machine-repairman  system.  The  graphics  display  indi¬ 
cates  the  model's  activities.  As  shown  in  Figure  7,  the 


1 .  Setting  the  Random-Number  Generator  Seed 

To  set  the  seed,  press  the  <s>  key.  The  program 
will  then  respond  by  asking  for  the  value  of  the  seed.  The 
seed  is  entered  by  typing  any  integer  in  the  range  of  1  to 
2,147,483,646  (=  231  -  2)  and  pressing  <ENTER>.  The  pro¬ 

gram  will  accept  a  non-integer  value;  however,  during 
processing  the  number  will  be  truncated  to  an  integer,  with 
the  result  possibly  being  different  than  the  desired  value. 
After  the  seed  has  been  entered,  the  program  returns  to  the 
program  menu. 

It  is  not  necessary  to  set  the  random-number  gener¬ 
ator  seed  to  run  the  simulation.  There  is  a  default  seed 
embedded  in  the  generator's  program  code,  and  it  is  updated 
automatically  with  each  call  for  a  random  number.  If  the 
simulation  is  stopped  and  then  restarted  without  setting 
the  seed,  the  generator  will  use  as  the  initial  seed  the 
value  remaining  from  the  last  call  in  the  last  simulation 


2 .  Changing  the  Model  Parameters 


When  <c>  is  pressed,  MACHREPR  will  successively  ask 
for  the  desired  values  of  the  mean  time  to  failure,  mean 
time  for  repair,  and  the  number  of  repairmen.  The  quanti¬ 
ties  are  entered  by  typing  a  positive  integer  and  pressing 
<  ENTER  >  in  response  to  each  question.  For  best  results, 
the  values  for  the  mean  time  to  failure  and  the  mean  time 
for  repair  should  be  at  least  ten.  The  program  will  accept 


parameters  that  can  be  altered  by  the  user  are  the 
helicopter's  mean  time  to  failure,  the  repairmen's  mean 
time  for  repair,  and  the  number  of  repairmen  to  employ. 
Additionally,  the  initial  seed  for  the  random-number 
generator  can  be  set. 


Figure  6:  The  Program  Menu. 


The  simulation  starts  and  ends  with  the  program  menu. 
Whenever  the  user  stops  the  simulation  or  whenever  the  sim¬ 
ulation  clock  reaches  the  preset  limit  of  9950  time  units, 
the  program  returns  to  the  menu.  The  simulation  then  can 
be  restarted  with  perhaps  a  difference  set  of  model  parame¬ 
ters,  or  the  program  can  be  ended  altogether. 
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Figure  5:  The  Title  Screen. 

3 .  Title  Screen 

The  display  will  remain  dark  for  a  few  seconds 
while  the  program  initializes  variables,  arrays,  and  func¬ 
tion  definitions.  When  initialization  is  done,  MACHREPR 
will  present  the  title  screen  as  shown  in  Figure  5. 


E.  THE  PROGRAM  MENU 

The  program  menu,  depicted  in  Figure  6,  appears  next 
and  provides  the  means  of  controlling  the  simulation.  With 
the  options  available  on  the  menu,  the  user  may  choose  to 
read  the  program  instructions,  change  the  model  parameters 
or  use  the  default  model  parameters.  The  program  instruc¬ 
tions  that  are  available  on  the  screen  are  a  summarized 
version  of  those  found  in  this  manual.  The  model 


command  to  make  an  identical  copy  of  the  distribution 
diskette . 

2  .  Starting  the  Simulation 

If  the  first  option  was  used,  insert  the  bootable 
diskette  into  drive  A  and  turn  on  the  microcomputer.  After 
a  few  seconds,  the  DOS  prompt  A>  will  appear  on  the  screen. 
It  is  here  that  the  batch  file  MACHREPR.BAT  can  be  used. 
Simply  type  'machrepr'  <ENTER>  to  start  the  simulation. 

If  using  the  non-bootable  diskette,  then  turn  on 
the  microcomputer  with  the  DOS  diskette  in  drive  A.  Next 
load  Advanced  BASIC  into  memory  by  typing  'basica'  <ENTER> . 
After  the  computer  is  in  BASIC's  command  mode,  remove  the 
DOS  diskette  from  drive  A  and  insert  the  diskette  contain¬ 
ing  the  simulation.  Type  'load  "machrepr . bas ", r '  <ENTER> 
to  bring  the  program  into  memory  and  to  start  the 
simulation . 

If  desired,  the  simulation  can  be  started  from  a 
second  disk-drive.  Begin  as  before  by  turning  on  the 
microcomputer  with  the  DOS  diskette  in  drive  A  and  loading 
BASICA.  With  the  application  diskette  inserted  into  drive 
B,  type  'load  "b :machrepr . bas "  ,  r  '  < ENTER > .  During  execu¬ 

tion,  MACHREPR  reads  the  program  diskette  and  expects  to 
find  it  in  the  default  disk-drive.  When  the  program  fails 
to  find  the  required  files  on  the  DOS  diskette  in  drive  A, 
it  asks  the  operator  for  the  correct  drive  label.  Simply 
respond  to  the  program's  question  by  pressing  <b> . 
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diskette  should  be  stored  in  a  safe  place  and  used  only  to 
make  a  replacement  for  an  unserviceable  application 
diskette.  Although  there  are  several  ways  that  the  user  may 
choose  to  create  an  application  diskette,  this  manual  will 
suggest  just  two: 

*  A  self -starting  (bootable)  diskette  containing  the 
operating  system  files,  the  BASICA  command  file,  and’ 
the  simulation's  program  files. 

*  A  simple  backup  copy  which  contains  only  the  program 
files  and  which  must  be  used  in  conjunction  with  a 
system  diskette. 

Owners  of  microcomputers  with  a  single  disk-drive  may  find 
the  first  option  more  convenient. 

To  make  a  bootable  diskette,  follow  the  instruc¬ 
tions  in  the  PC-DOS  manual  and  format  a  blank  diskette 
using  the  FORMAT  command  with  the  /S  parameter.  The  /S 
parameter  causes  the  computer  to  transfer  the  operating 
system  files  from  the  DOS  diskette  to  the  newly  formatted 
diskette.  The  Advanced  BASIC  file,  called  BASICA.COM,  is 
also  found  on  the  DOS  diskette  and  can  be  copied  onto  the 
application  diskette  with  the  COPY  command.  Finally,  use 
the  COPY  command  again  to  transfer  the  simulation's  six 
program  files  from  the  distribution  diskette  onto  the 
bootable  diskette. 

If  the  second  option  is  desired,  then  format  a 
blank  diskette  without  the  /S  parameter.  This  procedure 
does  not  transfer  the  operating  system  files  to  the  appli¬ 
cation  diskette.  The  next  step  is  to  use  the  DISKCOPY 


39 


APPENDIX  B 


MACHREPR . BAS  FLOWCHART 
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Repair  Module 
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Repair  Module  -  Par t  A 


Repafr  Module  -  Part  3 


Reoatr  Module  -  Part  C 
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SI, AX 

0060 

88 

FA 

MOV 

01,01 

006F 

01 

EO 

SAL 

AX,  1 

0071 

01 

02 

RCL 

01,1 

0073 

01 

EO 

SAL 

AX,  1 

0075 

01 

02 

RCL 

DX.l 

0077 

01 

EO 

SAL 

AX,  1 

0079 

01 

02 

RCL 

OX,  1 

007B 

2B  C6 

SUB 

AX, SI 

0070 

IB 

07 

SBB 

OX, 01 

007F 

E9 

012F  R 

JMP 

UNIF01 

; done  since  AIX<M 

0082 

33 

OB 

I 

HULT48:  XOR 

BX,BX 

;«ultiply  X  by  7A5 

0084 

88  CB 

MOV 

CM! 

0086 

01 

EO 

SAL 

AX,  1 

;by  nultiplymg  X  by  8 

0088 

01 

02 

RCL 

01,1 

0080 

01 

03 

RCL 

Bl,l 
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APPENDIX  D 


RNGEN.SRT  SOURCE-CODE  LISTING 

;  RNGEN.SRT  -  BASIC  USR  Routine 
;  by  B.Q.  Shubert 
;  Modified  by  R.E.  Neisen 
;  2  Hay  1984 
i 

;  This  subroutine  generates  and  returns  a  unitors 
;  randos-nusber  in  BASIC's  single-precision  forsat. 
;  The  seed  T  is  not  imported  nor  returned  to  BASIC, 
;  but  it  is  autosaticaily  updated.  The  algontha 
;  used  is: 

;  X  =  AX  sod  N 


;  with  A  =  7A5  and  H  =  2A31  -  1.  RNGEN  is  called 
;  froa  BASIC  using  the  stateients: 

i  RNGEN  =  0 

i  DEF  SEG  =  &H1A00 

;  CALL  RNGEN  (U!) 

;  U!  is  the  returned  randoa  nuaber  and  sust  be 
;  defined  prior  to  the  first  call  to  RNGEN. 

0000  CSEG  SEGMENT 

ASSUME  CS: CSEG, OS: CSE6, £S: NOTHING 


0000 

ORG 

0 

0000 

RNGEN  PROC 

FAR 

0000 

55 

PUSH 

BP 

;save  for  BASIC 

0001 

8B 

EC 

NOV 

BP,  SP 

jpoint  to  args  on  stack 

0003 

IE 

PUSH 

OS 

;save  for  BASIC 

0004 

8C 

C8 

NOV 

AX ,  CS 

;iake  data  accessable 

0006 

8E 

D8 

NOV 

OS, AT 

0008 

A1 

0161  R 

! 

NOV 

AX , LOSEED 

; get  old  T 

0008 

8B 

16  0163  R 

NOV 

DX , HISEED 

000F 

8B 

FO 

NOV 

SI, AT 

0011 

8B 

FA 

NOV 

01 ,  DX 

0013 

83 

FA  01 

CNP 

DX,  1 

;check  if  AIT  <  N 

0016 

77 

6A 

JA 

NULT48 

; 48  bit  sultiply  if  not 

0018 

72 

05 

JB 

NULT32 

;else  32  bit  sultiply 

001A 

3D 

F31D 

CNP 

AX,0F31DH 

001D 

73 

63 

JNB 

NULT48 

001F 

D1 

EO 

NULT32:  SAL 

AX,  1 

uultiply  T  by  7A5 

0021 

D1 

02 

RCL 

DX,  1 

0023 

D1 

EO 

SAL 

AT,  1 

002* 

01 

02 

RCL 

DX,  1 

0027 

D1 

EO 

SAL 

AX,  1 

0029 

D1 

02 

RC! 

DT.l 
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20147  LOCATE  8,1:  PRINT  *  The  numbers  listed  under  the  column  heading  'Lint'  are  the  luiting  va 
lues  as  ti»e  goes  to  infinity  of  the  above  defined  quantities.  These  values  are  computed  with  t 
he  balance  eguations  for  the  aodel.' 

201 48  LOCATE  15,1: PRINT  *  The  values  listed  under  the  column  heading  'Estimate'  are  the  estimate 
s  of  the  liaiting  values.  These  guantities  are  calculated  with  data  collected  frot* 

20149  PRINT  "the  simulation.  The  estimates  are  recom-puted  upon  every  occurance  of  an  event.' 

20150  LOCATE  25,2: PRINT  'Press  any  key  to  return  to  »enu...';:GQSUB  10000: GOTO  120 

20151  ' 

29999  '=====BASICA  AND  COLOR/GRAPHICS  ADAPTER  CHECK========= 

30000  DEF  SEG  =0: IF  (PEEK (fcH410)  AND  4H30)  <>  4H30  THEN  DEF  SE6:GQTQ  30007 

30003  LOCATE  3,1:PR1NT  'Sorry...' 

30004  PRINT  'You  do  not  have  the  color /graphics  monitor  adapter!’ 

30005  PRINT  'This  simulation  uses  graphics  and  requires  that  adapter.* 

30006  DEF  SEG  :  END 

30007  ON  ERROR  GOTO  30008:  PLAY  'pl6':E0T0  30011 

30008  WIDTH  80: CLS: LOCATE  3,1 

30009  PRINT  'This  siaulation  uses  advanced  BASIC.' 

30010  PRINT  ‘Reload  this  program  after  using  the  coaaand  ‘BASICA'.'iEND 

30011  ON  ERROR  GOTO  O:  RETURN 

30012  • 

30019  ===~FlLE-LOADING  ERROR  TRAPPING  RQUTINE============= 

30020  IF  ERR  <>  53  THEN  ON  ERROR  GOTO  0 

30022  CLSiLOCATE  3,1:PRINT  ’On  which  drive  can  this  program  be  found?" :GDSUB  10000 
30024  DRVE$=A$*V:  RESUHE  82 
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20044  LOCATE  ll,2:PRINT  ‘randoa-nuaber  generator  which  foraed' 

20045  LOCATE  13,2: PRINT  ‘the  basis  for  RNGEN.SRT.* 

20046  LOCATE  25,2:PRINT  ‘Press  any  key  to  continue. ..*;:SQSUB  10000:RETURN 

20047  ' 

20099  '  =====JNSTRUCT10NS==~=========“=s="=— ==-“====« 

20100  SCREEN  0,l:C0LQR  14,1 ,1:MDTH  40-.CLS 

20110  LOCATE  1,4:PRINT  ‘Instructions  for  the  Prograa  Menu*; 

20111  LOCATE  3,2:PR1NT  ‘Press  <s>  to  set  the  randoa-nuaber  seed';:PRINT  ‘and  enter  a  positive  intege 
r  when  asked. The  progra*  returns  to  the  «enu.‘ 

20112  PRINT: PRINT  *  Press  <c>  to  change  aodel  paraeeters.  The  progra*  dill  ask  for  the  aachines  ■ 
ean  tiae  to  failure  (HTTF),  the  repair-'; 

20113  PRINT  ‘aen's  aean  tiae  for  repair  (NTFRI,  and  the  nuaber  of  repairaen  to  eaploy.  The* 

20114  PRINT  ‘HTTF  and  the  HTFR  should  be  positive  integers  that  are  greater  than  10.  The  nuaber 

of  repairaen  can  be  1  or  2.  The  siaulation  starts  after  the  nuaber  of  repairaen  has  been  entered. 

« 

20115  LOCATE  18,2:  PRINT  ’Press  <d>  to  use  default  paraaeters.‘;:PRINT  ‘The  siaulation  Mill  iaaediat 
ely  start  hi th  the  values:  HTTF  =  15,  HTFR  =  10,  nuaber  of  repairaen  =  1.‘ 

20116  LOCATE  23,2:  PRINT  ‘To  end  the  prograa,  press  <e>.‘ 

20117  LOCATE  25,2:  PRINT  ‘Press  any  key  to  continue. G0SUB  10000 

20120  CLS: LOCATE  1,13:  PRINT  ‘keyboard  Coaaands'; 

20121  LOCATE  4,2:  PRINT  ‘During  the  siaulation,  you  aay  use  the*; : PRINT  ‘folloaing  keyboard  coaaands 
,  • 

20122  LOCATE  7,5:  PRINT  ‘<p>  -  pause  the  siaulation.* 

20123  LOCATE  9,5:  PRINT  '<c>  -  continue  the  siaulation.* 

20124  LOCATE  11,5:  PRINT  ’<s>  -  stop  the  siaulation  and  return* 

20125  LOCATE  13,11:  PRINT  ‘to  the  prograa  aenu.‘ 

20126  LOCATE  25,2:PRINT  ‘Press  any  key  to  continue... *;:60SUB  10000 

20130  CLS: LOCATE  1,13:  PRINT  ‘Graphic  Display'; 

20131  LOCATE  3,2:PRINT  ‘The  variable  naaes  in  the  display  are‘; sPRINT  'defined  as:'; 

20132  LOCATE  6,1:  PRINT  *TIHE:  current  tiae  on  siaulation  dock.’; 

20133  LOCATE  8,1:  PRINT  *N(t>:  no.  of  aachines  doun  at  tiae  t.‘; 

20134  LOCATE  10,3:  PRINT  *AR:  average  arrival  rate  of  aachines*; 

20135  LOCATE  11,7:  PRINT  ‘into  the  repair  queue.'; 

20136  LOCATE  13,3:  PRINT  'Nq:  average  waiting  tiae  of  aachines'; 

20137  LOCATE  14,7:  PRINT  ‘in  the  repair  queue.'; 

20138  LOCATE  16,3:  PRINT  *Lq:  average  length  of  repair  queue.*; 

20139  LOCATE  18,4:  PRINT  ‘N:  average  doan  tiae  of  a  Machine.*; 

20140  LOCATE  20,4:  PRINT  *L:  average  no.  of  doan  aachines.*; 

20141  LOCATE  22,2:  PRINT  'Avl:  aachine  availability.*; 

20142  LOCATE  25,2:  PRINT  ‘Press  any  key  to  continue. ..‘;:60SUB  10000 

20143  CLS: LOCATE  1,3:  PRINT  *Pd:  probability  that  a  down  aachine'; 

20144  LOCATE  2,7:  PRINT  ‘aust  aait  for  repair.*; 

20145  LOCATE  4,3:  PRINT  *Ps:  probilility  that  an  up  aachine'; 

20146  LOCATE  5,7:  PRINT  'aust  wait  in  cold  standby.*; 


15030  GOSUB  15038 
15032  1=131 
15034  GOSUB  15038 
15034  GOTO  15042 
15038  Y 1 =¥2-4 

15040  CIRCLE  11, Yl) ,12, 1,0, PI 
15042  ORAN  ’Cl  BU12  U4  L22  R44* 

15044  CIRCLE  (1,Y2) ,14,l,7.5»PI/6,11.5*PI/6 
15046  ORAN  *C1  B02  L4  NM-2,-15  M-8,-6* 

15048  ORAN  *011  M-4,+5  04  NL1  NR1  U3  M+6,-3* 

15050  ORAN  *8R20  M+6,+3  03  NL1  NR1  U4  M-4,-5  Ull* 

15052  ORAN  'M-8,+6  NM+2,-15  LA* 

15054  PAINT  STEP  (0,41,1,1 

15056  ORAN  *C0  L4  06  R1  U6  R6  06  R1  U4  L4* 

15058  CIRCLE  STEP  (0,4) ,1,0:PSET  STEP  (0,01,0 
15060  RETURN 

15062  'Statistics  Routine 

15064  LOCATE  14,22:PRINT  *HTTF: ‘j : LOCATE  14,32:PRINT  ’HTFR:'; 

15066  LINE  (161, 1131-1319, 113) ,2 

15068  LOCATE  16,22:PRINT  *Var':LOCATE  16,26:PRINT  ’Esti iate *: LOCATE  16,35:PR1NT  *Lnit* 
15070  LOCATE  1B,22:PRINT  *AR* 

15072  LOCATE  19,22:PR1NT  *Nq* 

15074  LOCATE  20,22: PRINT  *Lq* 

15076  LOCATE  21,22:PRINT  *N* 

15078  LOCATE  22,22:PRINT  *L* 

15080  LOCATE  23,22:PRINT  *Avl* 

15082  LOCATE  24,4:PRINT  *A*;:LOCATE  24,10:PRINT  *B*; 

15084  LOCATE  24,17:PRINT  *C,;:LQCAT£  24,22:PR!NT  *Pd*i 

15086  LOCATE  25,22:PRINT  *Ps*; 

15087  GET  (0,0)  -  (319,199) ,HA1NSCRN 

15088  GET  (1, 1431-149, 1831, GRNHELO 

15089  RETURN 

15090  ' 

19999  -s==sjiTLE  SCREEH=======*=— ======================= 

20000  SCREEN  0,1:  COLOR  14,1,1:  NIDTH  40:  CLS 
20010  LOCATE  2,12:  PRINT  ’GRAPHIC  SIMULATION* 

20012  LOCATE  4,18:  PRINT  'of  the* 

20014  LOCATE  6,9  :  PRINT  'MACHINE  -  REPAIRMAN  MODEL* 

20016  LOCATE  9,14:  PRINT  "by  R.E.  Nelsen* 

20020  LOCATE  14,5:  PRINT  'Submitted  in  partial  fulfillment* 

20021  LOCATE  15,2:  PRINT  *of  the  requireients  lor  the  degree  of* 

20022  LOCATE  17,1:  PRINT  'Master  of  Science  in  Operations  Research* 

20023  LOCATE  19,4:  PRINT  *fro«  the  Naval  Postgraduate  School* 

20024  LOCATE  20,11:  PRINT  'Monterey,  California* 

20025  LOCATE  22,2:  PRINT  ’Advisors:  J.D.  Esary,  A.F.  Andrus* 

20030  LOCATE  25,2:  PRINT  ’Press  any  key  to  continue...*;: GOSUB  10000 

20040  CLS: LOCATE  1,14:  PRINT  *Acknouledq«ent* 

20041  LOCATE  5,6:PRINT  ’The  author  Mould  like  to  thank' 

20042  LOCATE  7 , 2: PRINT  'Associate  Professor  Bruno  Shuhert  for* 

20043  LOCATE  9,2:PRINT  ’the  donation  of  his  asseibly-language* 
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12116  FAINT  STEP  (0,41,2,2 

12117  DRAM  "CO  L4  D6  R1  U6  R6  D6  R1  U6  L4' 

12118  CIRCLE  STEP  (0,4) ,2,0:PSET  STEP  (0,0), 0 
12120  BET  (X, 143)-  STEP  (48,40) ,REDHELOj X1=X1  +  19 

12130  CIRCLE  (11,165) ,4, 1:PA1NT  (11,1651,3,1 

12131  DRAM  ‘Cl  BH-4,-2  NL3  H*7,+l  M-4,-2’ 

12132  LINE  STEP  (0.6)-  STEP  (2,1),3,BF 

12133  LINE  STEP  (1,0)-  STEP  (-4,6) , 1 ,8F 

12134  LINE  STEP  (1,0)-  STEP  (2,6) , 1 , BF 

12135  ORAN  *Cl  L3  8U?  L2  N-4,-1  U1  M*4,+l  R2’:PAINT  STEP  (-1,1), 1,1 

12140  GET  (1,143)-  STEP  (48,40) , REDHEL02: RETURN 

12141  • 

12209  =="=FAILURE  DRAN1NS  SUBROUT  1  NES====================== 

12210  PUT  (1,105) ,GRNHELO:PUT  (1,143) ,REDHELD: RETURN: '  No  repairman 

12211  ■ 

12220  PUT  (1,105! , GRNHELO: PUT  (1, 143) .REDHELO 

12221  FOR  1=1  TO  300:NEIT:PUT  (1,143) , REDHELO 

12222  PUT  11,143) , REDHELD2: RETURN:  '  Nith  repairean 

12223  ' 

12230  PUT  (1,143) , 6RNHEL0: PUT  (1,105) .BRNHELO: RETURN: '  Cold  standby  to  flying  status. 

12231  ' 

13009  =====REPAIR  NODULE  SUBROUT  I  NES============— ======== 

13010  PUT  (1,143) ,REDHEL02:PUT  (1,1431 ,6RNHELQ: RETURN: '  To  cold  standby 

13011  ' 

13020  PUT  (!,143),REDHEL02:  FUT  (I, 143) .6RNHELD 

13021  FOR  1=1  TO  200:NEIT:  PUT  <1,1431 .6RNHELQ 

13022  PUT  (1,105) , GRNHELO:  RETURN:  ‘Repaired  and  goes  flying. 

13023  ' 

13030  FOR  1=1  TO  200-.NE1T-.PUT  (1,143) .REDHELO:  PUT  <1,143),REDHEL02:  RETURN:  'Assign  a  repairean. 

13031  ' 

M999  =====NAIN  SCREEN  DRAN1NB  SUBROUTINE  =======«==«==== 

15000  LINE  (0,8)-(319,8),2:LlNE  (0, 101 -(319, 10) ,2 

15002  LINE  (31 ,23) -(31,80): FOR  V=27  TO  75  STEP  16:LINE  (27, Y) - (31 , Y) : NEXT 

15004  LINE  (31 ,80) - (281 ,80) 

15006  FOR  1=36  TO  281  STEP  5 

15008  IF  (1=81)  OR  (1=131)  OR  (1=181)  OR  (1=231)  OR  (1=281)  THEN  LINE  (X,B0)-(I,B5)  ELSE  LINE  11,80) 
-(1,83) 

15010  NE1T 

15012  LINE  (0, 100) -(319, 100 1 ,2: LINE  (0, 102)- (319, 102) ,2 

15014  LINE  (t6t,l02l-(161, 1991,2 

15016  LOCATE  1,8:PR!NT  "MACHINE  -  REPAIRMAN  MODEL’ 

15018  LOCATE  4,3:PRINT  ’3’:L0CATE  6,3:PRINT  ’2’:L0CATE  8,3:PRINT  *1’:L0CATE  10,3:PRINT  *0' 

15020  LOCATE  3,1:PR1NT  ’N(t) ’ -.LOCATE  12,32:PRINT  ’Tue:’; 

15022  Helocopter  drawing  routine 
15024  PI  =  3.  141593: 1=25: Y2=165 
15026  G0SU8  15038 
15028  1=78: Y2=  127 
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9998  '»«**§  SUBROUTINES 

9999  ' 

10000  DEF  SEG: POKE  iH6A,0:  'General  input  routine—™5—1 

10001  A*5  INKEYI:  IF  A*5"  THEN  10001  ELSE  RETURN 

10002  ' 

10010  DEF  SEG:  POKE  iH6A,0:  ‘Pause  routine———555—— 

10011  A*5  INKEYS:  IF  A$5V  OR  A$5'C*  THEN  RETURN  ELSE  GOTO  10011 

10012  ' 

10099  ^"“STATISTICS  PRINT  SUBROUTINE55555— 5555555555555 

10100  LOCATE  18, Y: PRINT  USING  FHTlIjLAHBOA! ; 

10101  LOCATE  19,Y:PR!NT  USING  FHTlIjNO!; 

10102  LOCATE  20,Y:PRINT  USING  FHTISjLQ!; 

10103  LOCATE  21,Y:PRINT  USING  FHT1$;NI; 

10104  LOCATE  22,Y:PRINT  USING  FHT1$;LI; 

10105  LOCATE  23,Y:PRINT  USING  FMT1*; AVL! ; 

10106  LOCATE  24,Y:PRINT  USING  FNTlSjPDONN!; 

10107  LOCATE  25,Y:PR1NT  USING  FHT 1 $ ; PSTDBY I ; 

10108  RETURN 

10109  ' 

10209  ‘55555 NEXT  EVENT  SUBROUTINES55555555555555555555555555 

10210  IF  NXTFAIL ! =FAILTINE! (II  THEN  FAILINGHELO$=  *A*  ELSE  IF  NXTFAIL ! =FAILT IHE ' (2)  THEN  FAILINGHELOS 
5,B"  ELSE  FAIL1NGHEL0*5‘C" 

10211  RETURN 

10220  IF  NXTREPR !=REPRTINE ! (1)  THEN  REPAIREDHELOS5'A"  ELSE  IF  NXTREPR ! =REPRTIf1E f (2)  THEN  REPAIREDHEL 
0S5*8'  ELSE  REPAIREDHELOJ=-C" 

10221  RETURN 

10222  ' 

U049  —TIC  HARK  SUBROUTINE5555555555555555555— 5555555 

11050  AZ=STATE+l i ON  AZ  GOTO  11051,11052,11053,11054 

11051  Y=75:LINE  (XX, Y)  -  STEP  (5,0) ,2:RETURN 

11052  Y=j9:LINE  (XX, Y)  -  STEP  (5,0) ,2:RETURN 

11053  Y=43: LINE  (XX, Y1  -  STEP  (5,0) ,2:RETURN 

11054  Y=27:LINE  (XX, Y)  -  STEP  (5,0) ,2:RETURN 

11055  ’ 

12099  '=====FIRST  FAILURE  SUBROUTINE555555555555555555555555 

12100  DEF  SEG  =*iHl AOO: CALL  RNGENIU! ) :A*=  FNEXPONT'  (U!  ,HU')+TIHE 

12102  IF  FAILINGHELQ$=*B*  THEN  X554:X1578:REPRTIHE!  {2)5A!:FAILTIRE!  (2) 5lE+3 1 -.GOTO  12106 
12104  IF  FAILINGHELOt5'C'  THEN  1=107: Xl=131 jREPRTIHE ! (3) =A ! :FAILTIHE ! (3I=IE+31 
12106  PUT  (X,  105) , GRNHELO 

12110  PI  5  3.141593 

12111  CIRCLE  (XI, 161), 12, 2,0, PI 

12112  ORAN  *C2  BU12  U4  N-12,*10  H-8,-4  DIO  BH+20,-16  K+12,+10  R6  N+4,+5* 

12113  CIRCLE  (Xl,165) ,14,2, 7. 5»PI 76,11. 5»P 1/6 

12114  DRAN  *C2  BD2  L4  NH-2,-15  fl-B,-6’:DRAN  'Dll  H-4,+5  D4  NL1  NR1  U3  H+6,-3* 

12115  DRAN  *BR20  H+6,+3  D3  NL1  NR1  U4  H-4,-5  Ull'.'DRAN  ^-8,4-6  NH+2,-15  L4‘ 


2300  REPRTIHE! (AX)  *  1E+31  :FAILTIME !  (AX)  *  1E+31 
2310  IF  NUMRHAN-2  THEN  NAITREPR4*  MID* t ' ABC* , A1 t 1) 

2315  STARTQUE(AX)  *  TIME 

2320  GOSUB  12210:  "Ground  flying  helo:  no  repairean. 

2330  RETURN 

2331  ’ 

2999  «=«REPAIR  SUBROUT  I «===============»================ 

3000  AX*  INSTRCABC*  ,REPAIREDHELO$) 

3010  X  =  FNHELOPOSIT (REPAIREDHELOI) 

3020  ON  STATE  GOTO  3100,3200,3300. 

3021  " 

3100  STDBYHELO*=  HID*CABC\AI,1) 

3110  FAILTIHE! (AX)  =  IE+31 
3120  REPRTIHE! (AX)  =  1E+31 

3130  60SUB  13010:  'Change  red  helo  to  cold  stand-by  helo. 

3 UO  STATE  *  STATE-1:  RETURN 
3141  " 

3200  DEF  SEG  *6H1A00:CALL  RNSEN(U') 

3210  FAILTIHE!  (AX)  *  FNEXPONT!  (11! , ALPHA!  CTIHE 
3212  REPRTIHE' (AX)  *  1E+31 

3220  GOSUB  13020:  "Change  red  helo  to  green  and  fly. 

3230  IF  NUHRHAN-2  THEN  STATE*STATE-1:  RETURN 
3240  AX*  INSTRCABC*, HAITREPR41 
3242  I  =  FNHELOPOSIT (NAITREPR4) 

3244  SUNNS  *  SUHNStTIHE-STARTSUE(AX) 

3250  GOSUB  13030:  "Change  red  helo  to  red  helo  */  repairean. 
3260  CALL  RNGEN(U!):REPRTIHE' (AX)*  FNEXPONT1 <U!,HU!)+TIHE 
3262  NAITREPRS*** 

3270  STATE=STATE-l:  RETURN 

3271  " 

3300  DEF  SEG  =VH1A00:CALL  RNGEN(U') 

3310  FAILTIHE! (AX)=  FNEXPONT1 (U! .ALPHA ! ) ♦TIflE 

3312  REPRTIHE! (AX)  =  1E+31 

3320  GOSUB  13020:  "Change  to  green  helo  and  fly. 

3330  A1X*  INSTRCABC*, NAITREPRI) 

3332  X*  FNHELOPOSIT (NAITREPRt):  "Helo  at  head  of  queue. 

3334  SL'HNB  =  SUHNfQ+T IME-STARTQUE (All) 

3340  60SUB  13030:  "Change  red  helo  to  red  helo  «/  repairean. 
3342  CALL  RNGENtU!) 

3344  REPRTIHE! (A1X)=  FNEXPONT! (U!,HU!)+TIHE 

3346  IF  NUHRHAN*2  THEN  NAITREPRI***:SOTfl  3390 

3350  ON  AX  60T0  3362,3372,3382  j"  Advance  the  queue. 

3360  "REPA1REDHEL0*  =  *A* 

3362  IF  NAITREPR***B*  THEN  NA1TR£PR$**C*  ELSE  NAITREPR***B* 

3364  GOTO  3390 

3370  "REPAIREDHELO*  =  *B* 

3372  IF  NAITREPR»=*A*  THEN  NAITREPR$**C*  ELSE  NAITREPR*=*A* 

3374  GOTO  3390 

3380  REPAIREDHEL04  *  *C* 

3382  IF  NAITREPR*=*A*  THEN  NAITREPRI=*B*  ELSE  NAITREPR*=*A* 
3390  STATE*STATE-1:  RETURN 


1302  LASTEVNT  =  NITEVNT 

1310  IE  MXTFAIL ! <=MXTREPR!  THEN  GOSUB  2000  ELSE  6DSUB  3000 

1311  ' 

1399  1 m  Update  statistical  estiaates  *« 

1400  LAMBDA f  =SUMFA IL/TIME 

1410  IF  NUHRHAN- 1  THEN  PDOMN ! = (SUMP ! (21+SUHP1 (31 ) /TIME: LQ* =PDOMN 1 +SUHP ! (3) /TIME  ELSE  LQ!=SUHP! (3)/TI 
HE:PD(WN'=LQ! 

1420  L'=(SUMP! (11+SUHP! (21+SUHP1 (21+SUMP! (3! +SUMP* (3)+SUMPI (3> )/TIHE 

1430  MQ' =SUM«Q/SUMFA1L:M(=MQ!+1/MU! : PSTDBY 1 =SUMP 1 (0) /TIME 

1440  AVL!=SUHP!  (0) / T IME+ (SUMP !  (D+SUMP!  (D+SUMP!  (2) ) / (TIHE+TIME+TIME) 

1450  Y=25:SQSUB  10100:  'Print  estiaates  to  screen 

1451  ' 

1490  A$=  INKEYt:IF  AM"  THEN  1500 

1491  IF  A$=’s’  OR  AM’S*  THEN  120 

1493  IF  AM’p*  OR  A$=’P*  THEN  GOSUB  10010 

1495  DEF  SEG:  POKE  4H6A,0:  'Clear  key  buffer 

1500  WEND  : 'No  stop  coaaand  so  go  back  to  line  1100. 

1510  CLS:  LOCATE  12,1:  PRINT  'Prograa  tiae  liait  reached.’ 

1512  PRINT  ’Press  any  key  to  return  to  aenu...’ 

1514  G0SU8  10000:  GOTO  120 

1515  ' 

1999  '—^FAILURE  M0DULE======================— ==s=«===== 

2000  SUHFAIL-SUHFAIL+ 1 : STATE  =  STATE  +  1 

2020  IF  FIRSTFAILURE  THEN  GOSUB  12100: FIRSTFAILURE=FALSE: GOTO  2130 
2030  AM  INSTR(*ABC’,FAILINGHEL041 
2032  1=  FNHELOPQSIT (FAILINGHELOI) 

2040  ON  STATE  GOTO  2100,2200,2300 

2041  ' 

2100  DEF  SEG  =*tHl AOO: CALL  RNGEN1U!) 

2110  REPRTIHE! (AZI  =  FNEIPONT! (U! ,HU! ) +TIME 
2112  FAILTIME! (AZ)=  1E+31 

2120  GOSUB  12220:  '6round  flying  held  4  assign  repairaan. 

2130  AZ-  INSTRI’ABC’.STDBYHELOI) 

2132  1=  FNHELOPQSIT (STDBYHELOi) 

2140  CALL  RNGENIU! ) 

2142  FAILTIME! (AZ)  =  FNEIPONT! (U! , ALPHA ! ) +TIHE 

2150  GOSUB  12230:  1  Hove  stand-by  helo  to  flying  status. 

2160  STDBYHELOM": RETURN 

2161  ' 

2200  IF  NUHRNAN  =2  THEN  2240 

2210  GOSUB  12210:  '6round  flying  helo:  no  repairaan. 

2220  NAITREPRJ  =  MI D* ( ’ ABC’ , AX ,  1 ) 

2222  REPRTIHE' (AZ)  =  1E+31 : FAILTIME !  (AZ)  =  1E+31 
2224  STARTBUE(AZ)  =  TIHE 
2230  RETURN 

2240  GOSUB  12220:  'Ground  flying  helo  4  assign  a  repairaan. 

2250  DEF  SEG  =4H1A00:  CALL  RNGEN (U ! ) 

2252  REPRTIHE1 (AZ)  =  FNEIPONT' (U! , HU ' l+TIHE 
2254  FAILTIHE!  (AZ)  =  1E+31:  NAITREPR$=" 

2260  RETURN 

2261  ' 
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520  LOCATE  12,37:PRINT  USING  FHT2$;TIHE; 

521  GOTO  1000 

522  ' 

899  ==PROGRAfl  END  ROUTINE^— ==*=============“========— 

900  SCREEN  0,0:W1DTH  80:COLOR  6:PRINT  'Prograa  ended.. .*:END 

998  ' 

999  '-—CLOCK  HOOULE===================================»== 

1000  FAILTIHE! (1)=1E+31:  'Calculate  1st  failure  tiae  (B  l  0 
1010  DEF  SEE  -4H1A00:CALL  RNGEN (U f ) 

1012  FAILTIHE!  (2)  =  FNEXPQNT!  (U! , ALPHA! )  ♦TINE 
1020  CALL  RNGEN(U! ) 

1022  FAILTIME! (3)=  FNEXPONT!  IU! .ALPHA! )+TIHE 
1030  REPRTIHE! < 1 )=lE+3 1 :  * 1E+31  5  'infinity* 

1032  REPRTIHE!  12ME+31 
1034  REPRTIHE! (3) =1E+31 
1036  STDBYHELOI-'A' 

1099  **  Deteraine  next  event  »* 

1100  WHILE  NXTEVNK9950:  'Sets  upper  liait  to  prga  run  tiae. 

1110  IF  NOT  STATED  THEN  1120 

1112  HXTFAIL !=  FNH1N! (FA1LTIHE! (II ,FAILTIHE! (2) ,FAILTIHE! (31) 

1114  NXTREPR!  5  1E+31 
1116  SQSUB  102 10 : GOTO  1100 
1120  IF  NOT  STATED  THEN  1130 
1122  NXTFAIL!5  IE+31 

1124  NXTREPR!5  FNH1N1 IREPRT1HE! (1) , REPRTIHE! (2) .REPRTIHE' 13)) 

1126  GOSUB  10220: GOTO  1180 

1130  NXTFAIL'5  FNHIN1 (FAILTIHE! (1) .FAILT1HE! (2) .FAILTIHE! (3) ) 

1132  NXTREPR!5  FNHIN' (REPRTIHE! (1) .REPRTIHE! (2) .REPRTIHE! (3) ) 

1134  GOSUB  10210 
1136  60SUB  10220 

1180  IF  NXTFAIL! <5NXTREPR'  THEN  1184 
1182  NXTEVNT  5  FNROUNDUP tNXTREPR ! ) : GOTO  1200 

1184  NXTEVNT  5  FNROUNDUP (NXTFAIL!) 

1185  1 

1199  '»»  Run  dock  and  update  tiae  chart  ft 

1200  DEF  SEG: POKE  IH6A.0:  'Clear  key  buffer 
1210  WHILE  T1HE  <  NXTEVNT 

1220  TIHE=TIKE>1 

1230  IF  T1H£<550  THEN  GOSUB  11050: XX=XX+5:F0R  I51  TO  900:NEXT: GOTO  1240 
1232  DEF  SEG  5  4H1A00 
1234  CALL  SCRNSHIFT 
1238  XX=276: GOSUB  11050 

1240  LOCATE  12,37:PRINT  USING  FMT2$; TIME; 

1241  AtMNKEYt  :IF  A$5"  THEN  1260 

1242  IF  A$="s‘  OR  A*5'S*  THEN  120 

1244  IF  AJ=*p*  OR  A$5'P*  THEN  GOSUB  10010 
1250  DEF  SEG: POKE  4H6A.0:  'Clear  key  buffer 

1260  WEND 

1261  ' 

1299  '#«  Update  display  and  deteraine  next  event  t» 

1300  SOUND  120,3: SUMP ' (STATE ) =SUHP ! (STATE) ♦T1HE-LASTEVNT 


159  DEF  SEG  =  iHlflOO 

160  A1  !=  INT  (AI/16777210I):A2!=  INT((AI-Al!H6777210l)/65536!) 

162  'Poke  the  seed's  upper  2  bytes  into  RNSEN's  seed  storage. 

163  POKE  IH164,A1':  POKE  IH163.A2! 

164  A#=AI-A  1 ! *  1 67772 1 0I-A2 ! *6 5536 ! : A 1 !  =  INTIAI/256) :A2!=  AI-A1H256 
163  ‘Poke  the  seed's  loiter  2  bytes  into  RNSEN's  seed  storage. 

166  POKE  IH162,A1 ! :  POKE  IH161,A2! 

167  GOTO  120:  ‘Return  to  aenu 

168  ' 

199  -=====PARA«ETERS  RGUTIHE==— =========================== 

200  COLOR  15,5,0:CLS:HSG*=*>>Ef»ter  a  positive  integer  only.<<* 

202  LOCATE  l,6:PRINT  '»  CHANGE  MODEL  PARAMETERS  »':PRINT  SSI 

203  LOCATE  4,1: INPUT  'Enter  Mean  Tiae  To  Failure. A! : IF  A!<=0  THEN  LOCATE  23,3:PRINT  MSG$;:LOCATE 
4,1: PRINT  STRINGS (40,32) ; : GOTO  203  ELSE  ALPHA'=1/A' 

204  LOCATE  6,1: INPUT  'Enter  Mean  Tiae  For  Repair. ..*,A!:IF  A'<=0  THEN  LOCATE  23,3:PR1NT  MSSI;:LOCATE 
6,1: PRINT  STRINGS (40,32);:60T0  204  ELSE  MU'=l/A' 

205  LOCATE  23,1:PRINT  STRINGK40,32);:HSSI=" 

206  LOCATE  8, 1: INPUT  'Enter  nuaber  of  repairaen  (1  or  2)...*,NUNRMAN:IF  (NUMRMANOl )  AND  (NUMRHAN02 
)  THEN  LOCATE  8,1:PRINT  STRINSI  140,32) j :GOTO  206  ELSE  250 

210  ALPHA* M/15: HU 1  =  1  / 10: NUMRMAN=  1 

250  NXTFAIL ' =0:NXTR£PR ! =0: FIRSTFAILURE=TRUE: NXTEVNT=0: TINE=0: SUMFAIL=0: LASTEVNT=0: STATE=0: XX=31 : SUNN 
8=0 

252  FOR  1=0  TO  3:SUHP! (I) =0: STARTQUE ( I ) =0: NEXT 
260  ' 

299  =====THEORETICAL  STATISTICS  R0UTINE=================== 

300  P! (1)=2»ALPHA!/HU! 

320  IF  NUMRMAN=1  THEN  P! (2>=P! (l)*P! (l):P! I3)=P! (2)*ALPHA!/HU!  ELSE  P! (2)=P! (1)«ALPHA!/HU!:P! (3)=P! ( 
21 »ALPHA! / IHU'+HU!) 

324  P!I0)=1/I1+P!(1)+P!(2)+P!C3)) 

325  P! (1)=P! (1)*P! (0) 

326  P! (2)=P! (2) *P ! (0) 

327  P! (3)=P! (3)*P! (0) 

330  IF  NUMRMAN=1  THEN  PD0«N'=P! (2)*P! (3) : LO ' =PDONN ' +P ' (3)  ELSE  LB'=P! (3) :PD0NN*=P! (3) 

340  L(=P! (1)+2»P! (2)+3*P! (3) 

350  LAMBDA !  =  (2*P!  10) +2*P!  1D+P!  (21HALPHA! 

360  HQ ' =LQ !  /LAMBDA ! 

370  W'=L! /LAMBDA* 

380  PSTDBY!=P*  10) 

390  AVL'=P!(0IMP'<l)+P'fl)+P!(2>)/3 
395  ' 

499  •xss-PRINT  MAIN  DISPLAY  ROUT  1  NE====================- 

500  SCREEN  1,0:C0L0R  9,0 

502  IF  FIRSTSCRN  THEN  GOSUB  15000: FIRSTSCRN=FALSE: GOTO  510 
505  PUT  (0,0) ,HA1NSCRN 

510  LOCATE  14 , 27: PR INT  USING  FHT2S;  1 /ALPHA!; 

511  LOCATE  14,37:PRINT  USING  FHT2S; 1/HU! ; 

512  Y=33:G0SU8  1010Q:  '  Print  statistics  to  screen. 
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APPENDIX  C 


MACHREPR . BAS  PROGRAM  LISTING 


0  'MACHREPR. BAS,  Ver  3.1,  R.  E.  Nelsen,  30  Aug  1984 

1  CLEAR  :  KEY  OFF:  CIS 

2  60SUB  30000:  'Check  for  BASICA  end  color/graphics  adaptor 

3  DEF I  NT  A-Z 

4  DIN  P! (3) , SUMP ! (3) ,GRNHEL0(26S) , REDHELO ( 2i8) , REDHEL02 ( 268 ) 

5  DIN  T1MECHRT ( 1609) ,HAINSCRN(8002) 

6  DIN  FAILTIME ! (3) , REPRT 1  ME ? (3) , STARTQUE (4) 

7  RN6EN=0 : TRUE=- 1 : FALSE=0: U ' =0: SCRNSH IFT=3B4 

8  FNT 1 1= 'III . If I * : FHT2$s*f III* : DRVEI='A: * 

9  F I RSTSCRN= TRUE : MAI TREPR*= " * : FA I L I NGHEL0$= " ■ : ST0BYHEL0*= " " 

10  REPAIREDHEL0I=":S6I  =  STRING* (40,223) 

50  DEF  FNEXPONT ! (A1 ! ,A2! 1  =  -  L0S(A1!1/A2! 

51  DEF  FNROUNDUP(A')  =  INT (A! )  +  1 

52  DEF  FNMIN * (A1 ,B! ,C!)  =  (-(A!<B!  AND  A!<C!)*A!)+I-(B'<A'  AND  B ! <C ! > *B ! >  +  (- (C ' <A!  AND  C'<B!)»C!> 

53  DEF  FNHELOPOSIT (AI)=-(  1NSTR ( “ABC*  ,AI)  =  t)  +  M  INSTR ( ’ ABC’ , A*>  =2) *54 )  +  ( - (  INSTR CABC’ , A$)  =3)  *107) 
55  ' 

60  GOSUB  20000:  ‘Title  screen 

80  ON  ERROR  GOTO  30020:  Load  randoi  nuaber  generator. 

02  DEF  SEG  =  SrHlAOO:  BLOAD  DRVE**’RNGEN. SRT’ ,0 

84  ON  ERROR  GOTO  0 

85  ' 

90  BLOAD  DRVE*+ “SCRNSHFT . SRT' , 384 :  ‘Load  screen  shift  subroutine 

91  * 

100  ‘=====MA1N  PROGRAM  MENU=======— ======— =============== 

120  SCREEN  0,1: COLOR  14,3.0:NIDTH  40:CLS 

121  LOCATE  7,14:PR1NT  ’PROGRAM  MENU’zPRINT  SGI 

122  LOCATE  9,1:PRINT  'Constructions.' 

123  LOCATE  11,1: PRINT  ’<C>hange  aodel  paraeeters,* 

124  PRINT: PRINT  '<0>efault  aodel  paraaeters.* 

125  PRINT: PRINT  '<S>et  the  randoa  nuaber  generator  seed.' 

126  PRINT: PRINT  *<E>nd  the  prograa.':  PRINT  SGI 

130  LOCATE  21,2:  PRINT  'Enter  your  selection...'; 

131  ’ 

140  GOSUB  10000: AZ-  INSTR ( ' ICDSEi cdse' , Al) 

141  IF  AZ=0  THEN  120  ELSE  ON  AZ  SOTO  20100,200,210,150,900,20100,200,210,150,900 

142  ‘ 

149  ‘==”=RAND0N  NUMBER  GENERATOR  SEED  ROUT INE============ 

150  COLOR  14,0 ,0:CLS 

151  LOCATE  1,3:  PRINT  '»  SET  RANDON  NUMBER  GENERATOR  SEED  *':  PRINT  SGI 

152  LOCATE  4 ,2: PRINT  'Peraissible  seed  values  are  integers' 

153  PRINT  '  in  the  range:  1  to  2147483646.' 

154  LOCATE  7,2: INPUT  'Enter  the  seed  value. ..*;AI 

155  IF  A»< 1  OR  At >21474836461  THEN  LOCATE  7,2:  PRINT  STRI NGI (39 , 32) :  GOTO  154 


OOEE 

Dl  D2 

RCL 

DX,  1 

OOFO 

Dl  D3 

RCL 

BX,  1 

OOF2 

2B  C4 

SUB 

AX, SI 

00F4 

IB  D7 

SBB 

DX ,  DI 

OOFo 

IB  D9 

SBB 

BX,CX 

OOFB 

BB  FO 

MOV 

SI, AX 

OOFA 

8B  FA 

MOV 

DI ,  DX 

OOFC 

SB  CB 

MOV 

CX,  BX 

OOFE 

D!  EO 

SAL 

AX,  1 

0100 

Dl  D2 

RCL 

DX,  1 

0102 

Dl  D3 

RCL 

BX,  1 

0104 

Dl  EO 

SAL 

AX,1 

0104 

Dl  D2 

RCL 

DX,  1 

0108 

Dl  D3 

RCL 

BX,  1 

010A 

Dl  EO 

SAL 

AX,  1 

010C 

Dl  D2 

RCL 

DX.) 

010E 

Dl  D3 

RCL 

BX,  1 

0110 

2B  C4 

SUB 

AX, SI 

0112 

IB  D7 

SBB 

DX ,  DI 

0114 

IB  D9 

SBB 

BX ,  CX 

0114 

Dl  E2 

SAL 

DX,  1 

0118 

Dl  D3 

RCL 

BX,  1 

01 1 A 

Dl  EA 

SHR 

DX,  1 

one 

03  C3 

ADD 

AX ,  BX 

0 1 1 E 

S3  02  00 

ADC 

DX,  0 

0121 

81  FA  7FFF 

CUP 

DXJFFFH 

0125 

77  03 

JA 

SUBTRM 

0127 

EB  04  90 

JMP 

UNIF01 

012A 

40 

! 

SUBTRM:  INC 

AX 

012B 

81  DA  8000 

SBB 

DX.8000H 

012F 

A3  0141  R 

? 

UNIF01:  MOV 

LOSEED.AX 

0132 

89  14  0143  R 

MOV 

HISEED,  DX 

0134 

Dl  EO 

SAL 

AX,  1 

0138 

DI  D2 

RCL 

DX,  1 

013A 

32  C9 

XOR 

CL, CL 

013C 

Dl  EO 

FNDEXP:  SAL 

AX,  1 

013E 

Dl  D2 

RCL 

DX,  1 

0140 

72  04 

JC  • 

EXPONT 

0142 

FE  Cl 

INC 

CL 

0144 

ED  F4 

JMP 

FNDEXP 

0144 

Dl  EA 

EXPONT:  SHR 

DX,  1 

0148 

Dl  D8 

RCR 

AX,  1 

014A 

83  SO 

MOV 

BL.80H 

014C 

2A  D9 

SUB 

BL.CL 

> 


; BX, DX, AX ,  now  contain  AIX 

;get  X=AJX/2A31  to  BX 

;and  X’=A*X(mod  2A31>  to  DX, AX 

calculate  X”  =  X’+K 

; see  if  X”  <  M 

jdone  if  so  since  X”  can’t  =  N 

;el se  subtract  II  by  adding  1 
;and  subtracting  2A3I 

; store  new  X 

; Di vide  X  by  2A31  to  make 
;a  umfors  (0,1). 

; Set  a  counter  to  zero. 

;Shift  left  and  count  the 
, "number  of  leading  zeros, 
lleading  1  found. 

;Not  found,  add  1  to  counter 
land  look  at  next  bit. 
xShift  back  to  correct  value 
land  suppress  leading  1. 
iCompute  normalized  exponent. 


BL,  DX,  and  AX  now  contain  a  single  precision,  uniform 
(0,1)  random  number  in  BASIC's  floating  point  format. 


014E 

IF 

014F 

SB  7E  06 

0152 

88  25 

0154 

88  55  01 

0157 

88  75  02 

015A 

88  5D  03 

015D 

5D 

015E 

CA  0002 

0161 

4130 

LOSEED 

0163 

00A8 

HISEED 

0165 

RNEEN 

POP 

OS 

jrestore  secuent  register 

MOV 

DI,[8P]+6 

; get  addr  or  variable  0 

MOV 

EDI  3 ,  AH 

; pass  LSB  of  lantissa 

MOV 

I DI 3+ l , DL 

; pass  USB  of  eantissa 

MOV 

EDI  3+2, DH 

; pass  USB  of  eantissa 

HOV 

EDI  3+3, BL 

;pass  exponent 

POP 

BP 

jrestore  BP  for  BASIC 

RET 

“i 

; FAR  return  to  BASIC 

DM 

4130H 

; storage  for  1 

DM 

OABH 

ENDP 

APPENDIX  E 


SCRNSHFT . SRT  SOURCE-CODE  LISTING 


jiKHiHiKHItfHitililiKHIHItttlfiHlili 

{SCRNSHFT. SRT,  Ver  2  --  BASIC  USR  Routine 
;by  R.E.  Nelsen  --  27  July  1984 
i 

; Thi s  subroutine  shifts  5  coluans  to  the  left 
; the  4  rows  of  pixels  corresponding  to  the 
;aodel  s  state  vs  tiae  graph.  It  is  called 
;froa  BASIC  using  the  coaaands: 


;  SCRNSHFT  = 

0 

;  DEF  SES  = 

IH1A00 

;  CALL  SCRNSHFT 

J 

;This  subroutines  calls  the  RON  BIOS  Type  10H 

{interrupt  (Video  I/O)  using  routines  12  A  13. 

0000 

CSEG 

SEGMENT 

ASSURE  CS: CSEG 

DS; NOTHING 

{Establish  a  file  header  for  use  by  the 

{BASIC  BLOAD  coaaand. 

0000 

HEADER; 

0000 

FD 

DB 

OFDH 

{Code  for  BLOAD  file 

0001 

0000 

DN 

0 

;Seg  addr  location 

0003 

0000 

ON 

0 

{Offset  location 

0005 

0088 

ON 

RTN.LEN 

{Routine  length 

0007 

5 

SCRN.SHFT  PROC 

FAR 

0007 

55 

PUSH 

BP 

{Save  for  BASIC 

0008 

FB 

STI 

{Enable  interrupts 

0009 

B9  0024 

NOV 

Cl, 38 

;Scrn  coluan  36 

oooc 

BA  0018 

A6AIN:  NOV 

01,27 

;Scrn  row  27 

000F 

B4  00 

NOV 

AH, 13 

{Read  pix  at  ton  DI 

0011 

CO  10 

INT 

10H 

;ir  col  Cl.  Pix  in  AL 

0013 

83  E9 

05 

SUB 

Cl, 5 

{Shift  left  5  col  s 

0016 

B4  OC 

NOV 

AH, 12 

{Set  =  12  to  write 

0018 

CD  10 

INT 

10H 

{pixel  to  new  pos'n 

001A 

83  Cl 

05 

AOD 

Cl, 5 

{Shift  right  5  col's 

001D 

83  C2 

10 

ADD 

01,16 

{Row  43 

0020 

B4  00 

NOV 

AH, 13 

{Set  =  13  to  read 

0022 

CD  10 

INT 

10H 

{Read  pixel 

0024 

83  E9 

05 

SUB 

Cl, 5 

{Shift  left  5 

0027 

B4  OC 

NOV 

AH, 12 

{Set  =  12 

0029 

CD  10 

INT 

I0H 

(Write  pixel 

0026 

83  Cl 

05 

ADD 

Cl, 5 

{Shift  right  5  col  s 

83 


; Now  repeat  procedure  tor  retaining  too  rows. 


002E 

83  C2  10 

ADD 

01,16 

; Rom  59 

0031 

B4  0D 

nov 

AH, 13 

0033 

CD  10 

IMT 

10H 

0035 

83  E9  05 

SUB 

CX  ,5 

0038 

B4  0C 

NOV 

AH, 12 

003A 

CD  10 

INT 

10H 

003C 

83  Cl  05 

ADD 

Cl, 5 

003F 

83  C2  10 

ADD 

D  X ,  1 6 

{Row  75 

0042 

84  OD 

NOV 

AH,  13 

0044 

CD  10 

INT 

10H 

0046 

83  E9  05 

SUB 

Cl, 5 

0049- 

B4  OC 

NOV 

AH, 12 

004B 

CD  10 

INT 

10K 

004D 

83  Ct  06 

ADD 

CX,6 

{Shift  to  new  coluen 

0050 

81  F9  OUA 

CNF 

CX , 282 

;Is  last  col  reached? 

0054 

72  B6 

JB 

AGAIN 

;If  no,  do  next  col. 

9 

{Finished  with  rows,  so  now  lust  blank  oat  last 

;5  columns  of  each  ro«. 


i 


0056 

41 

INC 

CX 

; Cal  283 

0057 

B4  OD 

NOV 

AH,  13 

{Set  to  read  dot 

0059 

CD  10 

INT 

10H 

{Read  blank  dot 

005B 

8A  D8 

NOV 

61,  Al 

{Save  color  code 

005D 

B9  0115 

NOV 

CX,277 

{Col  277 

0060 

BA  001B 

NEN.CQL:  NOV 

DX  ,27 

{Row  27 

0063 

B4  OC 

NOV 

AH,  12 

(Set  to  write  dot 

0065 

BA  C3 

NOV 

AL,  BL 

{Set  backgrnd  color 

0067 

CD  10 

INT 

10H 

;Nrite  blank  dot 

0069 

83  C2  10 

ADD 

DX ,  16 

;Do  again  for 

006C 

B4  OC 

NOV 

AH, 12 

;row  43, 

006E 

8A  C3 

NOV 

AL ,  BL 

0070 

CD  10 

INT 

I  OH 

0072 

83  C2  10 

ADD 

DX,16 

;and  row  59, 

0075 

B4  OC 

NOV 

AH, 12 

0077 

8A  C3 

NOV 

AL,BL 

0079 

CD  10 

INT 

10H 

007B 

83  C2  10 

ADD 

DX ,  16 

;and  row  75. 

007E 

B4  OC 

NOV 

AH, 12 

0080 

BA  C3 

NOV 

Al,BL 

0082 

CD  10 

INT 

10H 

00B4 

41 

INC 

CX 

{New  col  on  right 

0085 

81  F9  OUA 

CNP 

CX,282 

{Last  coium? 

0089 

72  D5 

JB 

NEN.COL 

;No,  so  do  next  col 

i 

{Done,  so  recover  BP  and  return  to  BASIC. 
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